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Abstract 

Both  size  and  cost  of  mechanical  and  optical  inertial  navigation  systems  (INS)  have 
prevented  their  use  in  many  applications.  By  developing  a  micromachined  gyroscope  and 
combining  it  with  existing  micromachined  accelerometer  designs  to  form  a  low  sensitivity 
INS,  the  cost  would  be  reduced  by  a  factor  of  10  or  more.  The  lower  per  unit  cost  would 
open  new  areas  for  INS  use.  A  second  advantage  of  the  micromachined  INS  would  be 
its  size.  A  micromachined  INS  would  weigh  only  a  few  ounces  and  take  significantly  less 
space.  This  would  allow  more  room  for  electronics  and  the  weight  reduction  would  lead  to 
longer  mission  times. 

In  this  thesis  three  micromachined  gyroscopes  were  developed  and  fabricated:  a 
micromachined  comb-drive  vibrating  gyroscope,  a  planar  gimballed  gyroscope,  and  a  vi¬ 
brating  ring  gyroscope.  The  micromachined  gyroscopes  were  fabricated  using  the  Multi 
User  MEMS  (Micro  Electro  Mechanical  Systems)  Process  (MUMPS).  These  devices  were 
then  tested  in  AFIT’s  microelectronics  laboratories. 
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DESIGN  AND  FABRICATION  OF  A 


MICROMECHANICAL  GYROSCOPE 

I.  Introduction 

In  order  to  calculate  the  position  of  a  moving  body  accurately,  it  is  necessary  for 
navigation  systems  to  measure  both  the  motion  and  orientation  of  the  body.  An  inertial 
navigation  system  (INS)  can  be  used  to  provide  the  necessary  information  to  calculate  po¬ 
sition.  An  INS  typically  consists  of  three  orthogonally  oriented  accelerometers  and  three 
gyroscopes,  also  orthogonally  oriented.  Accelerometer  outputs  are  used  to  calculate  ac¬ 
celeration,  while  gyroscope  outputs  are  used  to  calculate  angular  orientation  of  the  body. 
Errors  in  calculating  angular  orientation  cause  the  INS  calculated  position  to  drift  with 
time,  making  it  necessary  to  update  the  INS  with  the  correct  position  or  to  develop  a 
means  of  modeling  the  drift  rate  in  order  to  maintain  accurate  position  data.  The  im¬ 
pact  of  the  drift  rate  on  positional  accuracy  makes  it  necessary  to  have  highly  accurate 
gyroscopes.  Sensitivity  constraints  of  devices  currently  used  make  them  costly  and  thus 
unsuitable  for  many  applications.  It  has  been  shown  that  by  integrating  position  infor¬ 
mation  from  the  Global  Positioning  System  (GPS)  with  a  medium  accuracy  INS,  accurate 
long-term  position  calculations  can  be  made  [3-6].  The  GPS  system  provides  highly  accu¬ 
rate  long  term  position  information.  On  the  other  hand,  high  frequency  maneuvers,  poor 
satellite  coverage,  and  the  time  necessary  to  calculate  position  causes  the  GPS  system,  as 
currently  implemented,  difficulty  in  tracking  highly  dynamic  movements.  However,  when 
accompanied  by  an  INS  capable  of  calculating  accurate  short-term  positions,  the  resulting 


integrated  navigation  system  can  provide  not  only  accurate  short-term  position  calculations 
(INS)  but  also  accurate  long-term  positions  (GPS).  The  accurate  long-term  characteristics 
of  GPS  make  it  possible  for  an  inexpensive,  moderately  precise  (rather  than  high-precision 
high-cost)  INS  to  be  integrated  with  GPS  to  provide  an  INS  capable  of  aerospace  naviga¬ 
tion.  One  example  of  an  integrated  system  is  the  integration  of  the  Litton  LR-80  attitude 
heading  and  referencing  system  (AHARS)  with  GPS.  The  LR-80  is  a  10  nmi/hr  system  in 
the  0.5  degree  heading  category,  a  poor  INS.  However,  when  it  is  integrated  with  GPS,  the 
overall  system  accuracy  is  20  m  absolute  position  error  and  0.25  degree  heading  error  [3]. 

Both  size  and  cost  of  mechanical  and  optical  inertial  navigation  systems  have  pre¬ 
vented  their  use  in  many  applications.  By  developing  a  micro  machined  gyroscope  and 
combining  it  with  existing  micromachined  accelerometer  designs  to  form  a  low  sensitivity 
INS,  the  cost  would  be  reduced  by  a  factor  of  10  or  more.  The  lower  per  unit  cost  would 
open  new  areas  for  INS  use.  A  second  advantage  of  the  micromachined  INS  would  be  its 
size.  The  Litton  LN94,  currently  used  in  the  F15  navigation  system,  weighs  over  45  lbs 
and  is  mounted  in  a  14  in  x  14  in  x  7  in  case.  A  micromachined  INS  would  weigh  only  a  few 
ounces  with  a  size  closer  to  3  in  x  3  in  x  0.5  in.  The  size  reduction  would  allow  more  room 
for  electronics  and  the  weight  reduction  would  lead  to  longer  mission  times.  The  microma¬ 
chined  INS  may  also  be  integrated  with  the  work  currently  being  accomplished  at  the  Air 
Force  Institute  of  Technology  (AFIT)  to  build  a  completely  jamming/spoofing- resistant 
GPS  receiver  on  a  printed  circuit  board  (PCB).  The  resulting  integrated  system  would 
possess  both  the  size  and  cost  advantages  mentioned  above  and  maintain  the  accuracies 
required  for  many  aerospace  applications. 
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Micromachined  gyroscope  designs  have  been  constructed  using  vibrating  rings  [7],  vi¬ 
brating  combs  [8] ,  and  planar  vibratory  rate  gyroscopes  [9] .  These  designs  are  implemented 
using  costly  custom  fabrication  processes.  If  a  standard  fabrication  procedure  could  be  used 
to  implement  a  design,  the  per-unit  cost  would  decrease  dramatically.  This  thesis  research 
will  be  limited  to  the  development  of  a  micromachined  silicon-based  gyroscope  fabricated 
using  the  Multi  User  MEMS  (Micro  Electro  Mechanical  Systems)  Process  (MUMPS)  and 
sensing  circuitry  developed  using  the  Metal  Oxide  Semiconductor  Implementation  System 
(MOSIS)  process. 

Development  of  a  micromechanical  gyroscope  capable  of  being  integrated  into  an 
INS  would  reduce  the  size  and  cost  of  current  inertial  navigation  systems.  By  using  a 
commercial  polycrystalline  silicon  (polysilicon)  micromachining  process,  availability  would 
be  increased  and  cost  of  functional  devices  would  be  further  reduced.  The  goal  of  this 
research  is  to  develop,  fabricate  and  test  a  micro  mechanical  gyroscope  using  a  commercially 
available  fabrication  process,  capable  of  sensing  an  angular  rotation  of  Ideglsec  with  a 
range  of  ±500deg / sec. 

Three  designs  are  investigated;  a  micromachined  comb-drive  vibrating  gyroscope,  a 
planar  gimballed  gyroscope,  and  a  vibrating  ring  gyroscope.  Each  of  the  designs  is  fabri¬ 
cated  and  tested  for  functionality  and  sensitivity.  All  three  of  the  designs  are  implemented 
using  the  MUMPS  process.  Sense  circuitry  is  fabricated  using  the  MOSIS  process.  In  this 
thesis,  fabricated  devices  are  compared  to  the  Charles  Stark  Draper  Laboratories  tuning 
fork  gyroscope  [9]  and  the  University  of  Michigan  resonating  shell  gyroscope  [7,10].  Devices 
fabricated  were  first  tested  at  atmospheric  pressure  and  then  in  a  near  vacuum.  Resonant 
frequencies  were  established  using  the  network  analyzer  available  in  AEIT’s  microelec- 
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tronics  laboratory.  Device  layouts  for  surface  micromachined  MEMS  are  completed  using 
Cadence  software  and  are  fabricated  by  MCNC  corporation  of  North  Carolina  (MCNC). 
Circuit  layouts  were  completed  using  Magic  software  and  fabricated  using  the  MOSIS 
process. 
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II.  Background 


2. 1  Vibrating  Gyroscopes 

In  1851  a  French  scientist  named  Leon  Focault  studied  the  earth’s  rotation  through 
the  use  of  a  large  pendulum  [11].  Focault  noted  that  as  the  earth  rotated  the  plane  in 
which  the  pendulum  was  oscillating  processed  about  its  central  axis.  Focault  correctly 
reasoned  that  the  precession  of  the  plane  of  oscillation  was  caused  by  the  Coriolis  effect. 
The  following  example  is  intended  to  explain  the  Coriolis  effect  in  greater  detail. 


Figure  2.1  (a)  Declaration  of  the  e  and  p  frames,  (b)  location  of  the  pendulum  in  the  p 

frame. 

2.1.1  Example:  Coriolis  effect.  A  snapshot  of  a  swinging  pendulum  is  shown  in 
Figure  2.1(b).  The  position  vector  to  the  pendulum  is  defined  as  Rp  and  the  p  reference 
frame  is  rotating  at  a  constant  angular  velocity  d;  with  respect  to  the  e  frame.  Given  the 
definitions  below,  the  acceleration  of  the  pendulum  with  respect  to  the  e  frame  expressed 
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in  the  p  frame  ^Ap  will  be  derived. 


(2.1) 


0 
0 

The  velocity  of  point  P  with  respect  to  the  e  frame,  coordinatized  in  the  p  frame  is: 

-  - 

^Vp  =  —Rp  +  LOxRp  (2.3) 

The  acceleration  of  point  P  with  respect  to  the  e  frame,  coordinatized  in  the  p  frame  ^Ap, 
is: 
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Expanding  the  Coriolis  acceleration  term  yields: 


0 

Xp 

2ujx  Rp  — 2 

0 

X 

Vp 

ip 

2cj  X  Rp 


2ijj^Xp 


0 


(2.6) 


(2.7) 


Equation  (2.7)  shows  the  results  of  expanding  the  Coriolis  acceleration  term  for  this  exam¬ 
ple.  The  result  is  a  compound  acceleration  in  the  Xp  and  Yp  directions.  The  acceleration 
of  the  pendulum  mass  will  cause  the  pendulum’s  plane  of  oscillation  to  process  about  the 
central  axis,  in  this  case  the  Zp  axis.  Given  two  reference  frames  rotating  with  respect  to 
each  other,  when  a  particle’s  acceleration  is  calculated,  the  Coriolis  acceleration  will  cause 
a  coupling  between  two  vibrational  modes  in  the  structure.  In  this  example  the  pendulum, 
originally  oscillating  in  the  YpZp  plane,  begins  to  rotate  about  the  Zp  axis.  Vibrational 
gyroscopes  depend  on  this  coupling  between  vibrational  modes  in  order  to  measure  angular 
rates;  their  sensitivity  is  directly  related  to  the  magnitude  of  coupling  present. 

2.1.2  Types  of  Vibrating  Gyroscopes.  Some  structures  other  than  Focault’s  pen¬ 
dulum  illustrate  coupling  between  vibrational  modes  when  subjected  to  rotation.  Of  these 
structures,  vibrating  prismatic  beams,  tuning  forks,  dual  accelerometers,  planar  gyroscopes 
and  vibrating  shells  are  available  either  commercially  or  as  research  devices  [12].  Each  of 
these  devices  will  be  discussed,  noting  their  significant  strengths  and  weaknesses. 
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Figure  2.2  Prismatic  beam  gyroscope. 


2. 1.2.1  Vibrating  Prismatic  Beams:  Vibrating  prismatic  beams  are  similar 
to  the  Focault  pendulum  in  operation.  Prismatic  beams  use  coupling  between  two  identi¬ 
cal  vibrating  modes  to  sense  rotation.  Historically  these  devices  have  shown  sensitivity  to 
temperature  variations  and  spurious  vibrations  [13].  Figure  2.2  shows  one  possible  imple¬ 
mentation  of  a  vibrating  prismatic  beam.  Because  of  its  sensitivity  to  off-ajds  vibrations, 
the  vibrating  prismatic  beam  is  not  suited  for  use  in  an  INS. 

2. 1.2. 2  Tuning  Forks:  As  early  as  1958,  the  Sperry  Gyroscope  company 

had  developed  a  navigation  grade  tuning  fork  gyroscope  [14,15].  The  Sperry  gyrotron  was 
machined  from  a  single  block  of  aluminum  on  a  macroscopic  level.  The  tines  of  the  tuning 
fork,  shown  in  Figure  2.3,  are  resonated  in  anti-phase.  When  a  rotation  is  input  about  the 
central  axis,  a  sinusoidal  torque  is  generated  about  the  support.  This  causes  the  support 
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Input  rate  (co) 


Direction  of  tine  vibration 

< -  - > 


Figure  2.3  Tuning  fork  gyroscope. 

to  undergo  torsion.  By  measuring  the  torsion  of  the  support,  the  input  angular  rate  may 
be  calculated.  In  the  tuning  fork  gyroscope,  the  primary  vibrational  mode’s  resonant  fre¬ 
quency  is  not  the  same  as  the  secondary  mode’s  resonant  frequency;  this  limits  sensitivity. 
Another  major  difficulty  in  the  tuning  fork  gyroscope  is  that,  if  the  mass  centers  of  the  tines 
are  slightly  off  center,  a  torsion  identical  to  the  one  caused  by  the  Coriolis  effect  will  be  seen 
on  the  output.  In  spite  of  these  difficulties,  a  micromachined  tuning  fork  gyroscope  was 
successfully  developed  at  Carl  Stark  Draper  Laboratories  using  a  custom  fabrication  pro¬ 
cedure  [9].  The  gyroscope  design  was  implemented  using  two  micromachining  techniques. 
The  first  technique  used  single  crystal  silicon  electrostatically  bonded  to  glass.  The  second 
technique  used  electroformed  nickel.  The  single  crystal  silicon  was  bulk  micromachined 
using  a  reactive  ion  etch  (RIE),  and  boron  was  used  as  an  etch  stop.  The  resulting  gyro- 
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scope  demonstrated  sensitivities  of  5000  deg/hr  in  a  60  Hz  bandwidth,  and  0.19  deg/sec  in 
a  1  Hz  bandwidth.  The  electroformed  nickel  gyroscope  was  patterned  in  photo  resist  and 
plated  using  nickel.  Specific  sensitivity  data  was  not  provided  on  the  nickel  gyroscope.  It 
was,  however,  stated  the  single  crystal  silicon  implementation  performed  better  than  the 
nickel  implementation  [9].  These  devices  do  not  have  the  sensitivities  desired,  0.1  deg/sec 
in  a  1  Hz  bandwith,  for  avionics  applications. 


Input  rate  (©) 

Direction  of  Direction  of 


Coriolis  acceleration 


Figure  2.4  Dual  accelerometer  gyroscope. 

2. 1.2. 3  Dual  Accelerometers:  Dual  accelerometers,  illustrated  in  Figure  2.4, 
are  a  variation  of  the  tuning  fork  concept.  Two  accelerometers  are  placed  in  parallel  and 
vibrated  in  anti-phase.  When  the  device  is  subjected  to  rotation,  the  Coriolis  acceleration 
acts  perpendicular  to  the  plane  of  vibration  [7].  This  type  of  gyroscope  is  sensitive  to  mass 
imbalances  and  is  suseptable  to  spurious  vibrations.  The  dual  accelerometer  approach  has 
not  been  pursned  because  of  its  sensitivity  to  off-axis  vibrations. 

2. 1.2. 4  Planar  Gyroscopes:  A  planar  gyroscope,  illustrated  in  Figure  2.5,  is 
a  two-gimbaled  gyroscope.  The  outer  gimbal  is  subjected  to  a  continuous  oscillation.  When 
the  device  is  subjected  to  rotation,  about  the  sensitive  axis,  the  inner  gimbal  begins  to 
vibrate.  The  major  drawback  to  this  design  is  that  the  drive  and  sense  oscillations  are  often 
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Figure  2.5  Planar  gyroscope. 


not  at  resonant  frequencies,  making  it  less  sensitive.  Charles  Stark  Draper  Laboratories 
has  successfully  constructed  a  monolithic  silicon  planar  gyroscope  [16-18].  The  device  is 
sensitive  to  angular  rates  of  10  rad/sec  (573  deg/sec);  this  device  is  severely  limited  by  its 
resolution.  Device  construction  was  completed  by  using  a  bulk  micromachining  process 
and  then  electroforming  gold  on  top  of  the  device.  A  similar  gyroscope  design,  constructed 
completely  from  silicon,  formed  a  silicon  pyramid  in  the  center  of  a  thin  membrane  [19]. 
The  membrane  allowed  the  pyramid  to  be  electrostatically  driven  to  resonance  in  one 
direction.  An  input  rotation  caused  the  pyramid  oscillation  to  couple  into  the  identical 
secondary  resonant  mode.  To  date  the  maximum  rate  sensed  with  these  devices  is  300 
deg/sec.  If  a  perfectly  manufactured  device  is  fabricated  the  noise  equivalent  rotation  rate 
would  be  0.07  deg/sec.  This  device  currently  suffers  from  the  same  limitations  reported 
by  Charles  Stark  Draper  Laboratories. 


2. 1.2. 5  Vibrating  Shells:  In  1890  G.  H.  Bryan  found  that  thin- walled  res¬ 
onators  could  be  used  to  sense  angular  rotations  [20].  Bryan  discovered  that  the  vibratory 
modes  of  a  ringing  wine  glass  process  when  the  glass  is  rotated  about  its  stem.  Further¬ 
more,  Bryan  noted  that  the  precession  was  due  to  the  Coriolis  acceleration  and  that  the 


2-7 


Standing  wave 
Antinode  axis 
after  rotation 


Figure  2.6  Vibrating  shell  gyroscope  [11]. 


angular  gain  could  be  computed  based  on  geometry.  Delco  Systems  Operations  developed 
a  Hemispherical  Resonating  Gyro  (HRG)  that  operates  by  measuring  the  precession  of  the 
nodes  in  a  vibrating  shell.  The  HRG  was  machined  on  a  macroscopic  level  from  fused 
quartz  with  a  metallized  surface.  The  surface  was  then  highly  polished  to  reduce  mass 
imbalances.  Delco  claims  that  the  HRG  is  a  nearly  ideal  sensor  because  of  its  [11]: 


1.  Solid  state  reliability 

•  no  rotating  parts 

•  small  inertial  dissipation 

2.  Insensitivity  to  magnetic  fields 

3.  Nuclear  ride-through  capability 

4.  Negligible  acceleration  sensitivity 

•  Same  performance  in  the  earth’s  atmosphere,  as  in  space 

5.  No  temperature  control  necessary 
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The  University  of  Michigan  in  coordination  with  General  Motors  has  experimented  with 
a  vibrating  ring  gyroscope  constructed  by  electro-forming  nickel  on  top  of  a  silicon  wafer 
processed  with  a  standard  Metal  Oxide  Semiconductor  (MOS)  process  [7,10].  The  resulting 
gyroscope  demonstrated  sensitivities  of  0.5  deg/sec  in  a  10  Hz  bandwidth.  The  gyroscope 
was  electrostaticly  balanced  to  account  for  mass  imbalances  and  operated  at  self  resonance 
to  allow  for  drift  of  the  resonant  frequency  with  temperature.  The  electrostatic  balancing 
must  be  accomplished  on  each  device,  making  mass  production  difficult.  The  major  factor 
limiting  the  device  sensitivity  is  the  use  of  nickel  as  a  device  material.  The  nickel  device 
has  a  resonant  Q  of  2  X  10^,  in  comparison  with  polysilicon  devices  which  typically  have  a 
resonant  Q  on  the  order  of  1  X  10®  [7]. 

2.2  Micromachining 

Silicon  is  a  mechanical  material  with  diverse  properties  [2].  Silicon  has  a  similar 
modulus  of  elasticity  and  a  greater  yield  strength  than  stainless  steel,  but  a  density  less 
than  that  of  aluminum.  In  addition,  silicon  has  almost  no  mechanical  hysteresis,  making 
it  ideal  as  a  sensor  material.  Table  2.1  compares  silicon  to  other  materials  used  in  device 
construction. 


Two  methods  of  constructing  micro-electro-mechanical  systems  (MEMS)  with  a  sil¬ 
icon  substrate  are  bulk  and  surface  micro-machining.  The  following  discussion  is  an 
overview  of  these  two  methods  of  device  fabrication. 


2-9 


Table  2.1  The  Mechanical  Properties  of  Silicon  [2]. 


Material 

Knoop 

hardness 

mm^ 

Young’s 

modulus 

Density 

cm^ 

Thermal 

Conductivity 

w 

cm^C 

Diamond 

53 

7000 

10.35 

3.5 

20 

Silicon  carbide 

21 

2480 

7 

3.2 

3.5 

Silicon  nitride 

14 

3486 

3.85 

3.1 

0.19 

Iron 

12.6 

400 

1.96 

7.8 

0.803 

Silicon 

7 

850 

1.9 

2.3 

1.57 

Tungsten 

4 

485 

4.1 

19.3 

1.78 

Stainless  steel 

2.1 

660 

2 

7.9 

0.329 

Molybdenum 

2.1 

3.43 

10.3 

1.38 

Aluminum 

0.17 

0.7 

2.7 

2.36 

2.2.1  Bulk  Micromachining.  In  the  1950’s  researchers  discovered  that  certain 
chemicals  etch  at  faster  rates  along  different  crystal  planes  of  silicon  [1].  Chemicals  such  as 
potassium  hydroxide  (KOH),  ethylenediamine-pyrocatechol  (EDP),  and  cesium  hydroxide 
(CsOH)  etch  along  some  silicon  crystal  planes  significantly  faster  than  others.  Other 
chemicals  etch  at  equal  rates  in  all  directions.  Chemicals  that  etch  selectively  are  called 
anisotropic  etchants;  chemicals  that  etch  at  equal  rates  in  all  directions  are  called  isotropic 
(Fig.  2.7).  Anisotropic  etches  are  characterized  by  sharp  corners  that  may  be  vertical  or 


Silicon  dioxide 


Silicon  dioxide 


Silicon  substrate 


Isotropic  Etching  Anisotropic  Etching 


Figure  2.7  Isotropic  vs.  anisotropic  silicon  etching  [21]. 


angled  depending  on  crystal  plane  orientations.  Figure  2.8  shows  the  effects  of  etching  two 
different  silicon  wafer  surfaces  with  the  same  anisotropic  etchant.  Crystal  directions  are 
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shown  and  the  corresponding  Miller  index  is  provided  [22].  The  54.7"  angle  is  a  result  of 
the  orientation  of  the  <  111  >  crystal  plate  with  respect  to  the  <  100  >.  Similarly  the 
90"  angle  is  a  result  of  the  <  111  >  crystal  plane’s  orientation  with  respect  to  the  <  110  > 
surface.  Structures  developed  using  anisotropic  etchants  also  depend  on  the  shape  of  the 


J1 


<111> 


Silicon  Substrate 


Figure  2.8  Anisotropic  etching  of  two  different  silicon  wafer  surfaces  [l]. 


opening  in  the  mask  layer.  Table  2.2  lists  etch  results  for  differently  shaped  windows 
exposing  a  <  100  >  surface. 

Table  2.2  Window  Shape  and  Etch  Results  for  <  100  >  Silicon  Surface  [1]. 


Window  Shape 

Etch  Results 

Square 

Truncated  Pyramid 

Rectangle 

Rectangular  Pit 

Circle 

Pyramidal 

Arbitrary  poly 

Rectangular 

Anisotropic  etchants  may  be  chosen  based  not  only  on  the  crystal  direction  that 
etches  faster  but  also  on  the  etch  rate  of  different  materials.  For  example,  if  a  layer  of 
silicon  dioxide  {Si02)  were  grown  on  top  of  a  <  100  >  silicon  surface,  a  window  were  cut 
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in  the  Si02  with  the  shape  shown  in  Figure  2.9(a),  and  the  wafer  were  placed  in  a  KOH 
solution,  then  the  structure  shown  in  Figure  2.9(b)  would  result. 


Window  exposing 
<100>  crystal  plane 


(a) 

Silicon  Dioxide 


(b) 

Figure  2.9  (a)  Etch  window  shape  exposing  the  <  100  >  crystal  plane  as  seen  from  above, 

(b)  Resulting  structure  after  anisotropically  etching  the  structure  in  (a). 


Bulk  micro-machined  structures  may  also  be  defined  through  the  use  of  etch  stops. 
Etch  stops  are  used  to  slow  the  etch  rate  of  silicon  further  when  exposed  to  certain  chem¬ 
icals.  Two  common  types  of  etch  stops  are  boron  doping  and  electro-chemical.  Regions 
highly  doped  with  boron  (lx  10^°cm“®)  cause  the  etch  rates  to  slow  by  as  much  as  1,000 
times.  Electro-chemical  etch  stops  may  also  be  used.  One  way  to  implement  an  electro¬ 
chemical  etch  stop  is  illustrated  in  Figure  2.10.  The  wafer  to  be  etched  is  suspended  in 
a  KOH-water  solution.  The  PN  junction  shown  in  Figure  2.10  is  reverse  biased  as  long 
as  the  P  layer  remains  intact.  When  the  P  region  is  fully  etched,  current  begins  to  flow 
and  a  layer  of  oxide  is  formed,  stopping  the  etch.  This  technique  may  be  used  to  create 
diaphragms  2  )um  thick  [1]. 
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Anode 


Cathode 


Figure  2.10  Electro-chemical  etch  stop  technique  [1]. 


Through  combined  use  of  isotropic  and  anisotropic  etchants,  under-cutting  of  Si02 
layers,  and  etch  stops,  complex  structures  may  be  developed.  Examples  of  structures 
created  using  bulk  micro-machining  are  accelerometers,  micro  valves,  and  hot  plates  [2,23, 
24]. 

2.2.2  Surface  Micromachining.  In  surface  micromachining,  structures  are  devel¬ 
oped  on  top  of,  not  within,  the  substrate  (wafer).  A  spacer  or  sacrificial  layer  is  grown  on 
top  of  the  substrate,  then  a  mechanical  layer  is  grown  and  patterned  on  top  of  the  sacrificial 
layer.  Once  the  mechanical  layer  is  patterned,  the  sacrificial  layer  is  etched  away,  leaving 
a  suspended  structure.  This  process  of  etching  the  sacrificial  layer  is  called  releasing  the 
structure.  Figure  2.11  shows  an  example  of  a  cantilever  beam  created  using  surface  micro¬ 
machining.  Surface  micro-machining  is  not  limited  to  one  mechanical  layer  [25].  In  the 


2-13 


(e) 

Figure  2.11  (a)  silicon  wafer,  (b)  wafer  with  sacrificial  oxide  grown,  (c)  oxide  patterned, 

(d)  mechanical  layer  grown,  (e)  sacrificial  oxide  layer  removed. 


previous  example,  after  the  mechanical  layer  is  grown  and  patterned,  a  second  sacrificial 
layer  may  be  grown  and  patterned.  A  second  mechanical  layer  may  then  be  grown  and 
patterned.  Once  the  second  mechanical  layer  is  grown  and  patterned  the  whole  structure 
may  be  released.  By  increasing  the  number  of  mechanical  layers,  the  complexity  of  designs 
possible  increases.  The  cost  of  using  more  layers  is  an  increased  difficulty  of  fabrication. 
Surface  micro-machining  has  been  used  to  fabricate  many  devices  such  as  micro-motors, 
resonant  combs  and  vibrating  rings  [26,27]. 


2.2.3  MUMPS  Process  Overview.  The  Multi-User  MEMS  Process  (MUMPS)  is 
a  3-layer  polysilicon  process  with  a  minimum  design  rule  of  2  /rm.  Table  2.3  provides  an 
overview  of  the  layers  in  the  MUMPS  process.  Once  all  the  layers  have  been  deposited 
and  patterned  the  structure  must  be  released.  This  step  is  completed  by  placing  the  die 
into  a  hydrofluoric  acid  (HF)  bath  dissolving  both  sacrificial  oxide  layers  and  releasing  the 
devices.  Figure  2.12  shows  a  released  mechanical  structure  identifying  each  of  the  layers 
present. 


2-14 


Table  2.3  Overview  of  Films  in  the  MUMPS  Process  [28]. 


Layer  Name 

Film  Description 

Comments 

Nitride 

0.5  pm  silicon  nitride 

isolation  layer 

Poly  0 

0.5  pm  polysilicon 

ground  plane,  wiring 

1st  Oxide 

2.0  pm  oxide  (PSG) 

1st  sacrificial  layer,  0.75  pm  dimples 

Poly  1 

2.0  pm  polysilicon 

first  mechanical  layer 

2nd  Oxide 

0.5  /im  oxide  (PSG) 

2nd  sacrificial  layer 

Poly  2 

1.5  pm  polysilicon 

second  mechanical  layer 

Metal 

0.5  pm  evaporated  metal 

low  resistivity  contact,  reflective  (gold) 

Poly  2 
Poly  1 
Nitride 
Dimple 

Figure  2.12  Released  device  cross  section  [29]. 

2.3  Chapter  Summary 

This  chapter  presented  three  main  topics.  The  first  was  the  Coriolis  effect.  The 
Coriolis  effect  is  the  underlying  principal  that  all  vibrational  gyroscopes  rely  on,  it  accounts 
for  the  coupling  of  energy  from  one  vibrational  mode  to  another.  The  second  was  a 
discussion  of  different  types  of  vibrational  gyroscopes.  There  were  five  possible  designs; 
the  vibrating  prismatic  beam,  tuning  forks,  dual  accelerometers,  planar,  and  resonating 
shell  gyroscopes.  Where  applicable  the  performance  of  fabricated  devices  is  related.  The 
third  topic  was  micromachining.  An  introduction  to  both  bulk  and  surface  micromachining 
was  presented  along  with  a  specific  discussion  of  the  process  used  to  complete  this  thesis. 
In  the  following  chapter  specific  device  theory  will  be  derived  for  a  tuning  fork  gyroscope, 
a  planar  gyroscope,  and  a  resonant  shell  gyroscope. 
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III.  Theory 


Three  types  of  vibrational  gyroscopes  are  investigated;  this  chapter  describes  the 
theory  behind  the  operation  of  each.  The  three  devices  are  the  vibrating  comb  gyroscope, 
planar  or  dumbbell  gyroscope,  and  the  vibrating  ring  gyroscope. 

3.1  Vibrating  Comb  Gyroscope 

In  this  section  the  structure  shown  in  Figure  3.1  is  analyzed.  The  torque  equations 


Figure  3.1  Comb  drive  gyroscope  illustration. 


are  generated  and  differential  equations  for  output  oscillations  are  derived  for  a  single  input 
about  the  device’s  sensitive  axis.  A  complete  derivation  for  an  arbitrary  input  rotation  is 


3-1 


completed  in  Appendix  A.  The  derivation  will  be  completed  for  one  vibrating  member.  A 
vibrating  member,  shown  in  Figure  3.2,  consists  of  the  four  support  flexures,  plate,  and 
the  drive  and  sense  combs.  Because  of  device  symmetry,  the  derivation  for  the  second 
vibrating  member  is  identical  to  the  first.  In  this  analysis  the  masses  of  the  plate,  drive 


comb,  and  sense  comb  are  combined  and  denoted  as  Mp.  M  will  refer  to  the  mass  of  the 
four  flexures. 

Similar  to  both  the  Sperry  Gyroscope  [15]  and  the  Charles  Stark  Draper  Laboratories 
gyroscopes  [9],  the  proposed  design  requires  the  two  vibrating  members  to  be  vibrated 
in  anti-phase.  In  order  to  reduce  the  voltages  necessary  to  sustain  an  oscillation  the 
device  is  oscillated  at  its  resonant  frequency.  The  resonant  frequency  of  the  device  may  be 
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determined  using  Rayleigh’s  method  [27]: 


fr 


''sys 


2-k  L(Mp  +  0.3714M)J 


(3.1) 


where  Mp  and  M  represent  the  mass  of  the  plate,  including  drive  and  sense  fingers,  and 
the  mass  of  the  supporting  beams,  respectively. 

For  this  analysis  the  trusses  joining  the  flexure  segments  are  assumed  to  be  rigid. 
As  a  result,  the  equivalent  spring  constant  for  a  single  beam,  k,  is  broken  down  into  three 
parts;  ki,  the  spring  constant  due  to  the  residual  stress;  ^2,  the  spring  constant  due  to 
beam  bending;  and  k^,  the  spring  constant  due  to  beam  elongation.  Residual  stress  refers 
to  the  stress  internal  to  the  material,  it  is  a  function  of  growth  conditions,  kgyg,  the  system 
spring  constant,  is  the  sum  of  the  spring  constants  from  each  of  the  four  support  beams 
and  is  written  as  in  Equation  (3.2). 


ksys  —  4(^1  +  A;2  +  ks) 


(3.2) 


The  vibrating  members  are  assumed  to  have  extremely  small  lateral  displacements.  Under 
this  assumption  ks  is  several  orders  of  magnitude  smaller  than  ki  and  ^2-  As  a  result  k^  is 
assumed  to  be  negligible. 

The  spring  constant  due  to  residual  stress  ki  is  calculated  in  this  section.  Figure 
3.3  shows  an  illustration  of  the  the  variables  necessary  to  calculate  ki-  The  axial  force, 
P  in  Figure  3.3,  is  constant  with  deflection,  d  in  Figure  3.3,  because  elongation  will  be 
accounted  for  by  kz-  The  force  stress  relationship  is  given  in  Equation  (3.3)  and  the  force 
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Figure  3.3  ki  spring  constant  illustration. 


balance  as  related  to  geometry  is  shown  in  Equation  (3.4)  [30]. 


P  ~  ctrA  =  bhaR 


(3.3) 


P  L 


(3.4) 


where  b  represents  the  beam  width  and  h  is  the  beam  thickness  in  the  direction  of  the 
desired  motion.  aR  is  the  residual  stress  of  the  material  and  A  is  the  cross  sectional  area 
of  the  beam.  The  spring  constant  ki  is  defined  in  Equation  (3.5) 


k^  = 


(3.5) 


L 


Figure  3.4  ^2  spring  constant  illustration. 


By  taking  advantage  of  symmetry,  the  analysis  for  the  k-^  calculation  is  broken  into 
two  parts  (see  Figure  3.4).  For  a  rectangular  beam,  Equation  (3.6)  defines  the  displacement 
of  half  the  beam  [30]: 


''2''  ZEI 


(3.6) 


Where  E  is  Young’s  Modulus  and  the  polar  moment  of  inertia  for  a  rectangular  beam,  I, 
is  defined  in  Equation  (3.7) 


/  = 


hb^ 

12 


(3.7) 


Multiplying  this  equation  by  two  and  solving  for  ^  yields  Equation(3.8) 


k2  — 


12EI 


(3.8) 


Substituting  Equation  (3.7)  into  Equation  (3.8)  yields  the  final  equation  for  ^2  in  terms  of 
the  beam  width  (6),  length  (L),  thickness  (h),  and  Young’s  Modulus  (E): 


^2  — 


Ebh^ 


(3.9) 


Finally,  substituting  Equations  (3.9)  and  (3.5)  into  Equation  (3.2),  the  equivalent  spring 
constant  for  the  vibrating  member  is  given  by  Equation  (3.10): 

,  .bh  Ebh^  ,  . 

keys  =  (3-10) 

This  value  of  k^ys  is  used  in  Equation  (3.1)  to  calculate  the  vibrating  member’s  resonant 
frequency. 
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The  quality  factor  Q  of  this  structure,  assuming  that  most  of  the  damping  is  from 
air  flow  underneath  the  plate,  is  [27]: 

Q  =  +  0.3714M)A;,p,]=  (3.11) 

liAp 

Where  fi  is  the  viscous  damping  of  air  (1.81  X  10~^N  s  at  20"C  [31]),  Ap  is  the  surface 
area  of  the  plate,  and  dp  is  the  distance  between  the  plate  and  the  substrate. 

Each  of  the  vibrating  members  is  driven  by  the  center  electrode  in  Fig  3.1  with  a 
voltage: 

Vd{t)  =  Ep  +  Vdsin{ut)  (3-12) 

Where  Vp  and  are  scalar  magnitudes.  This  causes  the  vibrating  members  to  oscillate 
in  antiphase.  The  lateral  (x)  displacement  of  the  vibrating  member  is  found  by  applying 
Hooks  law,  Equation  (3.13): 


2/90 'I 

P  _  Qx  ) 

""  2 

^  ^  (If)  9vd{tf 

dt  2ksys  dt 


(3.13) 

(3.14) 

(3.15) 


In  Equation  (3.14),  C  is  the  drive  electrode  interface  capacitance  defined  by  Equation 


(3.16). 


C  = 


n  X  €bl(x) 

diQ 


(3.16) 
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where  n  is  the  number  of  drive  fingers  on  the  vibrating  member;  de  is  the  distance  between 
the  stationary  and  moving  electrode  (2fimy,  and  l{x)  is  the  length  of  overlap  as  a  function 
of  displacement  in  the  x  direction.  For  completely  lateral  motion,  is  a  constant,  in 
other  words  l{x)  is  a  linear  function  of  x.  When  the  device  experiences  an  input  rotation, 
this  partial  derivative  is  no  longer  a  constant  and  is  instead  a  nonlinear  function.  The 
transfer  function  between  the  phasor  drive  voltage  Vd  and  the  phasor  displacement  X  at 
resonance  may  be  written  in  terms  of  the  magnitude  of 


X 

dC 

Vd 

^  k 

^sys 

dx 

(3.17) 


This  information  may  be  used  to  design  the  resonant  characteristics  of  the  device  in 
Figure  3.1.  Once  this  is  done,  it  is  necessary  to  calculate  the  response  to  an  input  rotation 
about  the  device’s  sensitive  axis.  Figure  3.5  shows  a  simplified  view  of  the  device  in  Figure 
3.1  defining  the  basis,  system  origin,  position  vectors  to  the  vibrating  members  centers  of 
mass,  the  magnitudes  R  and  H  used  to  write  the  position  vectors  in  the  basis,  and  noting 
the  sensitive  axis. 

For  this  analysis  let  both  vibrating  members  be  subjected  to  a  sustained  oscillatory 
displacement  given  by  Equation  (3.18),  where  a;  =  27r/r  is  the  radian  resonant  frequency 
calculated  by  using  Equation  (3.1)  and  is  the  amplitude  of  the  oscillation. 


Ax  =  ilsin{u}t) 


(3.18) 
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I  Sensitive  Axis 


I 


Figure  3.5  Simplified  tuning  fork  gyroscope  diagram. 

where  Ax  is  the  change  in  lateral  position.  The  position  vectors  of  the  left  and  right 
vibrating  member’s  center  of  mass  may  then  be  written  as: 


Rcmi  =  ~{R  +  ^sin{ut))xg  -  Hzg 
Rcmr  — '  (^R  T  g  H Zg 


(3.19) 

(3.20) 


Where  Rcmi  and  Rcmr  represent  the  position  vectors  to  the  left  and  right  members’  centers 
of  mass,  respectively.  When  the  device  is  subjected  to  a  rotation  about  its  sensitive  axis, 
at  a  rate  a  compound  acceleration  is  applied  to  the  resonating  members.  Using  the 
same  approach  as  in  the  Coriolis  acceleration  example,  the  accelerations  of  the  left  and 


3-8 


right  member’s  mass  centers  are  found  to  be: 

Scmi  —  +  ul)sin{u)t)  +  ojlR]xg  —  [(jJoR  +  ojo^sin{LL)t)  +  2LVo9,U)cos{u)t)]yg  (3.21) 

a^rnr  =  [— —  ul)sin(ut)  —  (jjlR]xg  +  [oioR  +  ojo^sin{ut)  +  2Uo^lcos{ut)]yg  (3.22) 

The  acceleration  of  each  member  is  needed  to  calculate  the  torque  that  each  member 
applies  about  the  central  axis  of  the  device.  The  torque  equation  in  vector  form  is 

f=RxF  (3.23) 

F  is  defined  by  Newton’s  second  law  [32] 

F  =  ma  (3.24) 

m  is  defined  in  this  case  as  the  mass  of  the  vibrating  member  and  a  is  the  acceleration  of 

the  member.  Solving  the  torque  equation  for  the  left  and  right  vibrating  members  yields 
the  torque  due  to  each  member  about  the  sensitive  axis: 

Tcmi  =  Tn[[—H{ojo{R  +  ilsin[ut)  +  2Uoilucos{ojt)}]x  g 
+[FI{^lu)^sin{(jjt)  +  Oi>l[R  +  Usin{ujt)]}]yg 
+[cUo(i2  +  ilsin{ut)Y  +  2u>o{R  +  ^sin{ut)){D,u)cos{ujt)]zg]  (3.25) 
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Tcmr  =  'rn]\^H{u}o{R  +  ^sin{(jjt)  +  2uJo^uJcos(ut)}]xg 
-[H{flu)^sin{ut)  +  u>o[R  +  Q,sin{u;t)]}]yg 

+[wo(i^  +  0,sin{u)t)Y  +  2Uo{R  +  ^lsin{u}t)){Q,ucos{LJt)]zg]  (3.26) 

The  net  torque  about  the  central  axis  is  found  by  adding  the  torques  due  to  each  vibrating 
member.  The  resulting  torque  is  shown  in  Equation  (3.27): 

f  =  <l)Xg  +  ^tjg  +  2m[ujo{R  +  ^sin{u>t)f  +  2u3„{R  +  ^sin{u>t)){Q.Locos{u:t)\Zg  (3.27) 

This  equation  demonstrates  an  interesting  result:  the  net  torque  only  has  a  compo¬ 
nent  about  the  sensitive  axis.  There  is  no  component  about  either  of  the  other  two  axes. 
This  is  only  true  because  of  the  symmetry  of  the  device  about  the  central  axis.  If  the 
center  of  mass  of  either  member  were  a  different  distance  from  the  origin,  then  there  would 
be  components  of  torque  about  the  other  two  axis.  It  can  also  be  shown  that,  if  the  device 
is  symmetric  about  the  line  connecting  the  two  mass  centers,  then  the  torque  about  the 
center  axis  is  the  same  at  both  the  bottom  and  top  connection  points.  This  is  a  desirable 
result;  torque  of  different  magnitude  about  the  top  or  bottom  flexure  would  create  highly 
nonlinear  output. 

A  relationship  may  now  be  defined  relating  the  angle  A^,  defined  to  be  positive  in 
Figure  3.6,  and  the  torque  generated  about  the  sensitive  axis  given  in  Equation  (3.27). 
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This  relationship  is  defined  below  [14]: 

I{t)Ag  =  T  -  CAg  -  XAg  (3.28) 

where  I{t)  is  the  system’s  moment  of  inertia  about  the  central  axis,  C  is  the  damping 
coefficient  of  air,  and  A  is  the  torque  required  to  rotate  the  structure  through  one  radian. 
A  is  determined  by  solving  Equation  (3.29)  for  T  where  9  =  57.3'’  (1  rad),  J  is  the  polar 
moment  of  inertia  for  the  support  defined  in  Equation  (3.30)  [33],  is  the  length  of  the 
support,  and  G  is  the  sheer  modulus  for  polysilicon  defined  in  terms  of  Young’s  modulus 
(E)  and  Poissons  ratio  (z/)  by  Equation  (3.31). 


9  = 


TL, 

GJ 


J  =  hh^ 


16 

y 


3.36 


h 

b 


(3.29) 

(3.30) 


G  =  E{l-u) 


(3.31) 
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Equation  (3.28)  may  be  resolved  into  principal  axes  and  /(/)  may  then  be  replaced  by 
the  (3,3)  scalar  element  of  the  transformed  matrix  (/^).  As  a  result,  the  scalar  nonlinear 
time  varying  differential  equation  for  Ag  may  be  written  as: 

:■  C  ■  X  ,  T  ,  , 

+  Y^g  +  Y^^  ~  T  (3.32) 

■^z 

The  output  of  the  vibrating  comb  gyroscope  has  been  completely  defined  in  terms 
of  a  second  order  differential  equation.  In  addition,  the  resulting  equation  for  Ag  may  be 
directly  compared  to  other  rate  gyroscope  output  equations. 

3.2  Planar  Gyroscope 

In  this  section  the  planar  gyroscope,  shown  in  Figure  3.7,  is  analyzed.  The  planar 
gyroscope  operates  by  establishing  a  sustained  oscillation  about  the  device’s  outer  gimbal. 
When  the  device  is  subjected  to  an  input  rotation  a  output  oscillation  is  established  about 
the  inner  gimbal.  The  torque  equations  are  generated  and  differential  equations  for  output 
oscillations  are  derived  for  a  single  input  about  the  device’s  sensitive  axis.  A  complete 
derivation  for  an  arbitrary  input  rotation  is  completed  in  Appendix  B.  In  order  to  analyze 
the  torque  generated  about  the  output  axis,  the  torque  generated  by  each  of  the  four  mass 
centers  must  be  calculated.  The  four  position  vectors,  basis,  and  system  origin  are  defined 
in  Figure  3.8.  The  corresponding  position  vectors  are  given  in  Equations  (3.33-3.36). 
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Figure  3.8  Planar  gyroscope  block  diagram. 


= 


-D 

-R 


(3.36) 


The  gyroscope  is  subjected  to  a  sustained  oscillation  about  the  yg  axis.  This  oscillation  is 
defined  by  Equation  (3.37). 


Ay  =  ilsin{(x)t)yg 


(3.37) 
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where  Ay  is  the  change  in  position  about  the  y  axis.  For  an  input  rotation  of  about 
the  Zg  axis,  the  rotation  rate  of  each  of  the  mass  centers  is  given  in  Equation  (3.38): 


0 

0 


(3.38) 


This  results  in  a  rotation  rate  of  oj  with  respect  to  inertial  space  given  in  Equation  (3.39). 
This  rotation  rate  has  components  due  to  the  input  rotation  and  the  sustained  oscillation. 


0 

u}Q,cos{u)t) 


(3.39) 


Using  the  above  definitions,  the  accelerations  of  each  of  the  mass  centers  is  calculated  and 
given  by  Equations  (3.40-3.43): 


Ai  —  [Ruz  —  Du^Q,^cos^{ut)  +  Dul]xg  [-Du^  —  Rijjl]yg 

+[(— Dw  —  Rjuj^oj)ilcos{ojt)  +  Du}^sin{ujt)]zg  (3.40) 

A2  =  [-RuJz  -  D{u9,cos{ut)f  A  ^lD]ig  +  {-Du^  +  R<^l]yg 

A{{—DA  -f  Rbj,,u>)^cos{u)t)  -b  Duj^Q,sin{uit)\zg  (3-41) 

A3  =  [—RoJz  +  Du^Q,^cos^{u)t)  —  Dul]xg  -b  [Dd>z  +  Ri^^]yg 

+[(jDd;  -b  RuJz^)^cos{u)t)  —  Duj^Q,sin{ujt)]zg  (3.42) 
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(3.43) 


A4  —  [BiOz  +  D.u)^0,'^cos^{u}t)  —  DLJ^]xg  +  [Dt^z  —  R‘^z]yg 
-\-[{Du  —  Ru)^cj)Q.cos{Lot)  —  DLJ^Q,sin(ut)]zg 

The  torques  generated  by  each  of  the  mass  centers  are  given  in  Equations  (3.44-3.47): 

Ti  =  m[[D RuQ,cos{ut)  —  DRiv^flsin(Lot)  -f  R‘^u>^ojQ.cos[ut)]xg 
+[—D^uClcos{iJt)  -t-  D^tJ^^sin{Lot)  —  DRu)z^^cos(LJt)\yg 

4-[— D^d;^  -{■  R^Uz  —  DR{ujQ,cos{LJtyf']zg\  (3-44) 

T2  =  m[[— Z).Rwficos(wt)  +  DRbJ^ilsin{u!t)  -t-  R^u)z<^^<^os{u>t)]xg 
+[—D'^u^cos((jjt)  D'^(jJ^Q,$in(u>t)  +  DRLOzLjQ,cos[ujt)]yg 

+[-D^u}z  +  R'^Cjz  +  DR{Loilcos{u)t)y]zg]  (3.45) 

Ts  =  m[[D R6jQcos{ijjt)  —  DRju'^0,sin{ut)  -1-  R^LJzOJ^cos{ut)\xg 
CjQ.cos{ojt)  +  D^uj^ilsin{u)t)  —  DRu>z(^^cos{ut)]yg 

-\-[—D'^Uz  +  R^<^z  —  DR{uQ,cos{ut)f\zg\  (3.46) 

T4  =  m[[— D Ruji}cos{ojt)  -|-  DRu^ilsin{LJt)  -f  R^ujzU;^cos(ut)]xg 
+[—D^uQcos{ut)  -|-  D^iJ^Q,sin{(x>t)  +  DRu}zU;Q,cos{ujt)]yg 

-{-[-D^d)z  +  R^>^z  +  DR{uUcos{iX)t)f]zg]  (3-47) 


3-16 


The  net  torque  is  found  by  summing  each  of  the  individual  torques  at  the  origin.  The 
resulting  net  torque  is  shown  in  Equation  (3.48): 

Tnet  =  4m[[R'^oj^u}9,cos{ojt)]xg  +  [- D'^u)ilcos(LJt)  +  D'^uj^^lsin{cjt)]yg  +  +  R^u^]zg 

(3.48) 

It  may  be  seen  by  the  final  torque  equation  that  the  structure  as  shown  in  Figure 
3.7  is  not  the  best  structure.  The  net  torque  equation  shows  the  effect  of  placing  the  mass 
centers  off  a  line  connecting  the  outer  gimballs  results  in  a  reduction  in  the  net  torque.  A 
structure  with  a  single  dumbell  like  mass  placement  along  the  line  connecting  the  outer 
gimballs  of  the  device  would  result  in  a  larger  output  torque.  This  may  be  shown  by  letting 
D  in  Equation  (3.48)  equal  to  zero. 

3.3  Vibrating  Ring  Gyroscope 

The  theory  for  the  third  gyroscope,  the  ring  gyroscope,  was  first  developed  by  Michael 
Putty  in  his  dissertation  at  the  University  of  Michigan  [10].  The  devices  developed  by  Dr. 
Putty  were  constructed  by  electroforming  nickel  to  generate  high  aspect  ratios  and  allow 
for  plate  capacitors  to  be  used  to  drive  and  sense  the  flexural  modes  of  the  device.  In 
this  thesis,  the  study  is  limited  to  use  of  a  standard  polysilicon  surface  micromachining 
process  (MUMPS).  The  aspect  ratio  in  the  MUMPS  process  is  limited  requiring  new  drive 
and  sense  mechanisms.  The  basic  ring  structure  evaluated  in  this  thesis  is  the  same  as 
Dr.  Putty’s,  making  most  of  the  device  theory  the  same.  Changes  to  the  sense  and  drive 
mechanisms  make  the  calculation  of  the  electrical  potential  energy  unique.  As  a  result 
of  the  similarities.  Dr.  Putty’s  dissertation  is  used  extensively  as  a  reference.  The  ring 
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structure  to  be  analyzed  is  shown  in  Figure  3.9.  The  device  is  suspended  otf  the  substrate 


Base 

Ring 


Support  Spring 


Figure  3.9  Basic  structure  of  the  ring  gyroscope. 


by  the  central  support  post,  or  base,  and  then  driven  at  its  flexural  resonant  frequency. 
The  resonant  mode  of  interest  is  the  second  order  mode  resulting  in  two  elliptically  shaped 
resonant  patterns  shown  in  Figure  3.10.  In  order  to  describe  the  device  behavior  normal 
mode  equations  were  developed  by  applying  Lagrange’s  equation  (3.49): 


dt  \dqi)  dqi  ^  % 


(3.49) 


Where  T  is  the  kinetic  energy  of  the  system,  U  is  the  potential  energy,  and  Q  is  the 
generalized  force  found  from  the  amount  of  work,  dW,  to  perform  a  virtual  displacement, 
dq  [34].  The  analysis  results  in  a  set  of  coupled  second  order  differential  equations.  The 
analysis  relies  on  the  fact  that  a  complex  vibrational  displacement,  u,  may  be  written  as 
a  linear  combination  of  its  normal  mode  displacements,  with  complex  amplitudes, 
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(b) 

Figure  3.10  Flexural  vibrations  in  a  ring  gyroscope;  (a)  primary  vibrational  mode  (b) 
secondary  vibrational  mode. 

as  given  by  Equation  (3.50): 

n 

u  =  Y,^i^i  (3.50) 

i=l 

Vibrational  modes  other  than  the  two  fundamental  flexural  modes  contribute  only  minor 
second  order  effects;  as  a  result  they  are  not  included  in  the  device  model. 

It  is  now  necessary  to  define  the  coordinates  used  in  the  analysis:  Figure  3.11  defines 
the  X  and  y  directions;  Figure  3.12  defines  the  generalized  coordinates  that  give  the  complex 
amplitude  of  the  two  flexural  modes;  and  Figure  3.13  defines  the  generalized  coordinates 
used  to  calculate  the  generalized  force  Q*  in  Lagrange’s  equations.  The  first  generalized 
coordinate,  qi,  is  located  an  angle  of  6^  counter-clockwise  from  the  x  direction  and  the  q2 
direction  is  45  degrees  past  the  qi  direction.  The  angle  6t  defines  the  location  of  the  qi’ 
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direction.  The  vibrational  displacement  may  then  be  described  in  terms  of  its  x  and  y 


Figure  3.11  Definition  of  the  x  and  y  coordinate  directions. 


/ 


(a) 


0^+4^ 


Figure  3.12  (a)  definition  of  the  coordinate  direction  (b)  definition  of  the  q2  coordinate 

direction 


displacements,  and  Uy,  its  mode  shape  functions  for  the  x  displacement,  $3,1  and  $3.2, 
its  mode  shape  functions  for  the  y  displacement,  $j,i  and  $j,2)  and  coordinatized  in  terms 
of  the  generalized  direction  vectors  qi  and  52-  Equations  (3.51)(3.52)  define  the  and  Uy 
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(a)  (b) 

Figure  3.13  (a)  definition  of  the  coordinate  direction  (b)  definition  of  the  ^2’  coordi¬ 

nate  direction. 

displacements: 


Ux  — 

(3.51) 

«y  =  +  ^y2q2 

(3.52) 

To  facilitate  the  analysis  further,  the  ring  will  be  broken  up  into  nine  parameterized  seg¬ 
ments,  one  for  the  ring  itself  and  eight  for  the  support  springs.  Figure  3.14  shows  the 
parameterization  for  each  of  the  elements  in  terms  of  the  x  and  y  directions.  Applying 
these  parameter  definitions  to  Equations  (3.53)(3.54)  and  adding  the  parameter  p,  the 
final  X  and  y  displacement  functions  may  be  written  as  in  Equations  (3.53  and  3.54). 


u^{p)  ~  +  ^x2ip)q2  (3.53) 

Uy{p)  =  ^ylip)qi  +  ^y2{p)q2  (3.54) 
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•etememt location  x,(!d),y/.Q) 


y-axis  I  dlsplacfincot  tt,(8),n  (0) 


Figure  3.14  Ring  gyroscope  parameterization  [10].  In  the  upper  portion  of  the  figure 
the  parameter  6  represents  the  location  of  a  point  on  the  ring.  In  the  lower 
portion  6  represents  the  location  of  a  point  along  one  of  the  semi-circular 
support  springs. 
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In  order  to  apply  Lagrange’s  equations  to  this  system,  it  is  necessary  to  calculate  the 
kinetic  energy,  potential  energy,  and  generalized  force.  The  following  sections  complete 
these  calculations. 

3.3.1  Kinetic  Energy.  The  derivation  of  the  kinetic  energy  of  the  system  was 
completed  by  Michael  Putty  in  his  dissertation  [10].  The  definitions  necessary  are  provided 
along  with  the  results  of  the  analysis.  Let  the  ring  structure  be  located  in  a  reference  frame 
OXYZ  at  ra  =  XaX  +  yay  (see  Figure  3.15).  Let  the  OXYZ  frame  be  rotating  with  a  rate  of 
oj  =  (flj.,  fly ,  and  translating  with  respect  to  an  inertial  frame  OoXoYoZo  with  a  velocity 
of  Vo  =  Vxo&  +  VyoV  +  v^oZ-  The  kinetic  energy  is  determined  by  calculating  the  inertial 


Figure  3.15  Coordinate  definitions  for  kinetic  energy  derivation  [10]. 


velocity  and  finding  its  magnitude  squared,  multiplying  by  the  density  of  the  material  p. 
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and  integrating  over  the  volume  of  the  device.  The  resulting  equation  is  shown  in  Equation 


(3.55): 


(3.55) 


After  the  integration  is  performed,  the  equation  is  simplified  and  the  definitions  of  the 
small  X  and  y  displacements  are  substituted  into  the  equation.  The  result  is  given  in 
Equation  (3.56): 


T  —  -MiqI  +  -M^ql  +  ~  929i)  + 

+  2(^i^2  +  (3.56) 


where 


Mi  =  lp[^li  +  ^lx]dV 

M2  =  lp[^l2  +  ^l2]dV 

7  =  /  -  ^r2^y2]dV 

JV 

^1  =  j^P^lidV 

132  =  Ip^UdV 

a,=  f  p^l,dV 

Jv 

0^2=  I  P^l2dV 
Jv 

Pi  =  lp[^li  +  ^l]dv 
P2  =  lp[^l2+^l2]dV 
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It  is  important  to  note  that,  in  the  second  term  of  the  kinetic  energy  equation,  the  Coriolis 
coupling  depends  only  on  the  rotation  about  the  device  sensitive  axis.  It  should  also  be 
noted  that  the  kinetic  energy  is  independent  of  the  location  of  the  ring  and  the  translational 
velocity  of  the  moving  frame.  This  shows  the  independence  of  mounting  location  and,  to 
first  order,  insensitivity  to  linear  accelerations.  This  completes  the  kinetic  energy  analysis; 
it  is  necessary  next  to  derive  the  potential  energy  of  the  system. 

3.3.2  Potential  Energy.  The  potential  energy  terms  may  be  broken  up  into 
two  categories,  the  mechanical  energy,  11^-,  of  the  ring’s  differential  elements,  and  the 
electrical  potential  energy,  Ug.  The  mechanical  potential  energy  of  the  ring  gyroscopes 
designed  for  this  thesis  is  identical  to  the  mechanical  potential  energy  derived  by  Putty  in 
his  dissertation.  The  electrical  potential  energy,  however,  is  unique  to  these  devices.  The 
mechanical  potential  of  the  ring  is  given  by  Equation  (3.57). 

Urn  =  \Kxql  +  \K,ql 

”  j-  — FPwi — “ " 

K2  =  EI  [  ~  +  yp^x2  - 

Jp  [ip  +  J/p]  j 

The  electrical  potential  energy,  Ug,  is  analyzed  via  the  specific  electrode  sets  affecting 
the  ring.  Specifically,  the  electrical  potential  energy  may  be  separated  into  the  energy  from 
the  drive  electrodes  {U^),  sense  electrodes  (Hj),  and  from  the  remaining  electrode  sets  (Ug). 
The  electrical  potential  energy  is  dependent  on  the  size  and  shape  of  the  electrodes.  In 
the  case  of  the  ring  gyroscope  developed  at  the  University  of  Michigan,  the  ring  had  flat 


(3.57) 

(3.58) 

(3.59) 
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faces  interfacing  with  flat  electrodes.  This  was  possible  because  of  the  thickness  of  their 
structures,  25/rm,  compared  to  the  maximum  thickness  from  the  MUMPS  process  of  3.5^m. 
As  a  result,  the  drive  and  sense  electrodes  are  constructed  from  a  series  of  interdigitated 
combs  (Figures  3.16  and  3.17).  There  are  one  fixed  comb  and  a  comb  attached  to  the  ring 
that  make  up  a  single  electrode. 

The  potential  energy  of  the  system  may  be  expanded  and  written  as  a  linear  combi¬ 
nation  of  the  potentials  from  each  electrode  set;  this  is  shown  in  Equation  (3.60): 

Ug  =  Ud  Us  Ug  (3.60) 

Ud,  Ug,  and  Ug  are  replaced  with  their  definitions  in  equation  3.61 

c,  =  -  E  \v;c„  -  £  lv?C„  -  Y.  (3.61) 

n  =  l,5  ^  n=2,6  ^  n=3,4,7,8  ^ 

' - V - '  ' - V - ^  ' - V - ' 

drive  sense  ground 

In  Equation  (3.61)  the  subscripts  on  the  summations  refer  to  the  electrode  numbers  as¬ 
signed  in  Figure  3.16,  U„  is  the  total  voltage  applied  to  the  electrode,  and  C„  represents 
the  capacitance  generated  by  the  interaction  between  the  electrode  and  the  ring.  The 
derivation  of  each  capacitance  for  each  of  the  eight  electrodes  is  the  same,  as  a  result  the 
analysis  wiU  be  completed  for  one  electrode  and  then  substituted  into  Equation  (3.60). 
The  capacitance  is  calculated  by  assuming  infinitely  small  parallel  plate  capacitors  in  par¬ 
allel.  This  reduces  the  capacitance  equation  to  an  integral  form.  Figure  3.17  shows  a  single 
electrode  and  ring  interface.  Let  the  initial  finger  overlap  length  be  defined  as  L  and  the 
nominal  position  be  defined  as  dg.  The  change  in  capacitance  may  be  written  in  terms  of 
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3 


7 


Figure  3.16  Electrode  location  definitions. 

the  small  displacements  and  Uy.  Figure  3.18  shows  the  interface  between  one  ring  finger 
and  two  of  the  electrode  fingers.  It  is  necessary  to  integrate  along  the  finger  from  its  start 
to  the  end  of  the  overlap.  Let  the  end  of  the  stationary  electrode  finger  be  defined  as  zero. 
The  change  in  overlap  distance  is  then  defined  by  Equation  (3.62): 

\uj;\cos{d)  +  |uj,|co5(90‘’  —  9)  (3.62) 

the  magnitudes  of  the  x  displacement  and  the  y  displacement  must  be  found  by  applying 
the  law  of  cosines  because  of  the  non-orthogonality  of  the  qi  and  q2  directions.  In  general 
the  magnitude  of  may  be  written  as  in  Equation  (3.63)  and  the  magnitude  of  Uy  may 
be  written  as  in  Equation  (3.64): 
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Figure  3.17  Single  ring  gyroscope  electrode. 


Figure  3.18  A  single  ring  linger  and  its  stationary  electrode  counterparts. 


3-28 


1«.|  =  -  2$.i(0)$.2(^)co5(135°)  (3.63) 

\uy\  =  +  ^y2{ey  -  2^y,{e)^y2{0)cos{u5o)  (3.64) 

The  total  overlap  at  any  time  instant  then  becomes: 

L  —  \ux\cos{0)  +  |My|cos(90®  -  6)  (3.65) 


In  order  for  the  overlap  distance  to  be  used,  L  must  he  defined  in  terms  of  the  generalized 
coordinates  and  .  By  applying  the  law  of  sines  to  a  vector  of  magnitude  L  located  at 
an  angle  of  6,  the  vector  L  expressed  in  terms  of  the  generalized  coordinates  is  given  by 


Equation  (3.66): 


L  = 


Lsin{45°  +  0o  —  0) 
sm(135o) 


01  + 


Lsin{0  —  0o) 
5m(135'’) 


(3.66) 


and  the  capacitance  from  one  ring  finger  is  defined  by  the  integral: 


C  = 


Jo  .do  T  ^d 


dl 


(3.67) 


where,  e  is  the  permeativity  of  the  medium  and  t  is  the  layer  thickness,  do  is  the  initial 
stationary  to  free  electrode  spacing,  and  Ad  is  the  change  in  spacing  from  the  stationary 
to  the  free  elctrode.  The  first  term  inside  the  integral  comes  from  the  stationary  electrode 
above  the  moving  electrode  and  the  second  comes  from  the  lower  electrode.  If  either  of 
the  stationary  combs  is  missing  (e.g.,  at  the  ends  of  the  electrode),  the  corresponding  term 
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may  be  omitted.  When  both  stationary  fingers  are  present,  the  integral  simplifies  to: 


C  = 


2etd„ 

Jo  dl  +  lS.d? 


(3.68) 


A  McLaurin  series  expansion  is  used  about  Ad  to  linearize  this  equation.  This  simplifies 
the  integral  to: 


C  = 


(3.69) 


Evaluating  the  integral  yields  the  capacitance  due  to  one  moving  finger  surrounded  by  two 
stationary  fingers: 


C  = 


2etL 

~T 


(3.70) 


It  is  also  necessary  to  calculate  the  result  for  the  upper  stationary  electrode  being  present 
only,  and  for  just  the  lower  stationary  electrode  as  well.  The  result  of  completing  the 
analysis  for  the  upper  electrode  only  is: 


fi  fi  ft 

C^  =  -  dl--  Addl+-  /  Ad^dl  (3.71) 

d  Jo  d?  Jo  d^  Jo 

The  result  for  only  the  lower  electrode  being  present  is: 

ff  rt  ft 

C,  =  -  dl+—  Addl  -  —  /  AdJdl  (3.72) 

d  Jo  d?  Jo  d^  Jo 
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Where  Ad  is  defined  as  in  the  direction  of  increasing  0  and  is  related  to  the  and  Uy 
displacements  by  Equation  (3.73): 

|Ad|  =  |Mj;|sm(0)+  \uy\cos{9)  (3.73) 

In  order  to  maintain  consistency,  it  is  necessary  to  write  a  vector  Ad  in  terms  of  qi  and 
52-  This  is  given  in  Equation  (3.74) 

Ad  =  —  |Ad|cos(90°  —  d  +  9o)qi  +  |Ad|cos(45°  6  —  0o)q2  (^-^'i) 

The  capacitance  for  an  electrode  is  found  by  summing  the  individual  effects  of  the  ring 
fingers.  The  effects  are  dependent  on  time  and  the  angular  location  of  the  finger.  Using 
the  definitions  above,  the  capacitance  for  the  nth  electrode  may  be  written  as  shown  in 
Equation  (3.75). 

Cn  =  +  C,  (3.75) 

The  potential  energy  may  now  be  written  as  a  function  of  the  ring’s  mode  shapes.  The 
mode  shapes  for  the  analysis  must  come  from  a  finite  element  analysis  tool.  The  use  of 
these  tools  is  beyond  the  scope  of  this  thesis;  as  a  result  a  closed  form  solution  for  the 
potential  energy  of  the  system  is  not  included.  Investigation  of  a  closed  form  solution 
for  these  mode  shapes  and  a  finite  element  analysis  of  the  other  gyroscopes  investigated 
provides  an  excellent  opportunity  for  future  research.  The  Nastran  finite  element  package 
available  at  AFIT  would  be  an  ideal  package  to  use.  The  last  piece  needed  to  implement 
Lagrange’s  equations  is  the  generalized  force  Q. 
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3.3.3  Generalized  Force.  The  non-conservative  forces  rise  mainly  from  damping, 
both  material  damping  and  viscous  damping.  The  generalized  force  Q  is  a  function  of  the 
amount  of  work  dW  necessary  to  create  a  virtual  displacement  dq.,  and  is  defined  in  terms 
of  the  qi  and  ^2’  coordinates  by  Equations  (3.76)  and  (3.77): 


Qi'  = 

Q2'  = 


dW 

Qi' 

dW 

h' 


—Ciqif 

-C^q-il 


(3.76) 

(3.77) 


Cl  and  C2  are  the  independent  damping  coefficients.  In  order  to  maintain  consistency,  the 
qi  and  0^2 ’  coordinates  must  be  translated  into  the  qi  and  ^2  coordinates.  The  results  of 
this  coordinate  transformation  are  given  below: 

Qi  =  -[CiCOs^(2(/))  -b  C2sin'^{2(f))]qi 

-[(Cl  —  C2)cos{24>)sin{24>)]q2  (3.78) 

Q2  =  -[(Cl  -  C2)cos{24')sin{2<f>)]qi 

-[CiCos^{2^)  +  C2sin^{2^)]q2  (3.79) 


Where  ^  is  defined  as  : 


(l>  —  0T  -\-  do 


(3.80) 


The  damping  accounts  for  the  loss  of  energy  in  the  system  due  to  damping  both  in  the 
material  and  due  to  viscous  damping  in  the  air.  This  completes  the  evaluation  of  the  terms 
necessary  to  implement  Lagrange’s  equations. 
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In  this  section  the  ring  gyroscope  was  analyzed  using  a  work  energy  approach.  Be¬ 
cause  of  design  choices,  the  electrode  capacitance  is  not  a  simple  function  of  radial  dis¬ 
placement;  instead  it  is  a  function  of  a  linear  combination  of  the  mode  shape  functions. 
The  complexity  of  the  structure  makes  it  necessary  to  use  a  finite  analysis  program  to 
attain  the  mode  shape  functions.  The  difficulty  in  learning  these  tools  makes  this  analysis 
beyond  the  scope  of  this  thesis. 

3.4  Chapter  Summary 

In  this  chapter  the  theory  of  operation  for  the  three  gyroscopes  investigated  was 
developed.  Versions  of  each  gyroscope  were  designed  and  fabricated.  In  the  following 
chapter  each  of  the  designs  is  discussed  and  experimental  procedures  are  recorded. 
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IV.  Experimental  Approach  and  Procedures 

In  the  previous  chapter  a  theoretical  derivation  was  completed  for  the  operation  of 
three  types  of  mechanical  gyroscopes.  Each  of  the  gyroscopes  was  then  fabricated.  The 
purpose  of  this  chapter  is  to  describe  the  resulting  designs  and  their  testing.  The  designs 
will  be  discussed  in  the  order  of  creation.  There  were  four  fabrication  runs;  MUMPS  6, 
MUMPS  7,  MUMPS  8  and  MUMPS  9.  The  submission  date  of  each  run  is  given  in 
Table  4.1  below.  Sections  4.1  -  4.9  describe  the  designs  included  on  each  of  the  fabrication 

Table  4.1  Fabrication  Dates  For  MUMPS  Runs. 


Run  Name 

Due  Date 

MUMPS  6 

1  Dec  94 

MUMPS  7 

28  Feb  95 

MUMPS  8 

8  May  95 

MUMPS  9 

14  Aug  95 

runs  and  provide  the  statistics  specific  to  each  run.  The  variables  that  were  assumed 
constant  over  all  the  fabrication  runs  are  provided  in  Table  4.2.  The  device  descriptions 


Table  4.2  Common  Design  parameters. 


Parameter  Name 

Symbol 

Density 

P 

Young’s  Modulus 

E 

160  MPa 

Poisons  Ratio 

V 

0.23-0.3 

are  followed  by  the  experimental  set-ups  and  procedures  used  to  test  for  device  resonant 
frequency  characteristics. 

4.1  MUMPS  6 

The  MUMPS  6  fabrication  run  statistics  are  given  in  Table  4.3.  In  that  table,  the 
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Table  4.3  MUMPS  6  Fabrication  Statistics. 


Film 

Thickness  (A) 

Stress  (MPa) 

Resistivity  (ohm-cm) 

poly  0 

19.7  (C) 

1.84E-3 

poly  1 

5.1  (C) 

2.58E-3 

poly  2 

5.1  (C) 

2.72E-3 

metal  (Cr/An) 

30  (T) 

3.39E-6 

nitride 

6213 

19.2  (T) 

- 

oxide  1 

- 

- 

oxide  2 

- 

- 

stress  of  the  layer  is  followed  by  a  C  or  a  T;  C  represents  a  compressive  stress  and  T  a  tensile 
stress.  The  dashes  indicate  values  not  computed.  The  nitride  layer  acts  as  an  isolation  layer 
so  its  resistivity  is  very  high.  The  two  oxide  layers  are  the  sacrificial  layers  of  the  process. 
These  layers  provide  spacing  between  mechanical  layers  and  are  removed  once  the  structure 
is  released.  MUMPS  6  was  the  first  fabrication  run  used  in  this  thesis  and  device  theory 
was  under  initial  investigation.  There  were  two  classes  of  devices  created.  The  first  was  a 
vibrating  comb  gyroscope.  In  order  to  maximize  the  output  of  this  gyroscope,  a  multi-layer 
structure  was  investigated,  using  different  flexure  and  support  lengths.  Figure  4.1  shows 
the  device  as  created.  In  addition  the  legend  attached  to  all  Scanning  Electron  Microscope 
(SEM)  photos  is  labelled.  The  multilayer  structure  was  created  strictly  adhering  to  the 
MUMPS  process  design  rules  given  in  Table  2.3.  As  a  result,  the  device  is  not  symmetric 
in  the  vertical  direction.  A  cross  section  of  the  resulting  structure  is  shown  in  Figure  4.2. 
As  the  comb  gyroscope  theory  matured,  it  became  obvious  that  this  type  of  device  was  not 
ideal.  The  asymmetric  nature  of  the  vibrating  plates  lends  itself  to  twisting  under  normal 
operation.  As  a  result,  an  output  will  be  seen  for  a  zero  input  rotation.  In  addition,  the 
sense  capacitance  is  too  small  to  provide  a  reliable  signal.  There  were  also  initial  planar 
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Scale  Reference 


Figure  4.1  Initial  comb  gyroscope  design.  The  legend  attached  to  all  Scanning  Electron 
Microscope  Photographs  is  labelled. 


Gold  — 

Poly  2  - 
PlP2via 
Poly  1  - 


Silicon  Substrate 


Figure  4.2  Initial  comb  gyroscope  design  cross  section. 


gyroscope  designs  included  on  this  fabrication  run,  but  layout  errors  in  these  designs  made 
them  unusable. 


4.2  MUMPS  7 


The  characteristics  of  the  MUMPS  7  fabrication  run  are  given  in  Table  4.4.  Versions 


Table  4.4  MUMPS  7  Fabrication  Statistics. 


Film 

Thickness  (A) 

Stress  (MPa) 

Resistivity  (ohm-cm) 

poly  0 

4880 

- 

2E-3 

poly  1 

20100 

8.0  (C) 

1.65E-3 

poly  2 

14706 

11.1  (C) 

1.54E-3 

metal  (Cr/Au) 

5318 

21  (T) 

3.30E-6 

nitride 

6049 

41.0  (T) 

- 

oxide  1 

20400 

- 

- 

oxide  2 

5281 

- 

- 

of  the  tuning  fork,  ring,  and  planar  gyroscopes  were  fabricated  in  this  run.  Each  device  is 
discussed  in  terms  of  its  sense  mechanism,  drive  mechanism,  and  design  considerations. 


4.2.1  Vibrating  Comb  Gyroscope.  There  were  two  versions  of  the  vibrating  comb 
gyroscope  implemented  in  this  run.  The  first  is  named  combi.  This  design  uses  combs 
both  to  drive  the  sustained  oscillation  and  to  sense  the  output  oscillation.  Figure  4.3  shows 


4-4 


the  fabricated  structure.  Variations  on  this  device  consisted  of  changes  in  structural  layer 


Figure  4.3  Combi  gyroscope  design  in  MUMPS  7. 


and  alterations  in  the  width  of  the  central  support.  Table  4.5  lists  the  specific  variations 
and  design  geometry  information.  The  variations  in  the  dimensions  of  the  torsion  bar 
are  intended  to  allow  for  different  torsional  spring  constants  for  the  output  flexure.  The 
variations  should  result  in  different  device  sensitivities.  The  device  spring  constant  and 
resonant  frequency  are  calculated  by  applying  the  equations  developed  in  Chapter  3.  The 
only  device  variation  affecting  the  resonant  frequency  is  the  structural  layer.  The  resonant 


4-5 


Table  4.5  Variations  and  Geometry  Information  for  the  Comb  1  Gyroscope. 


Parameter 

Value(s) 

Structural  Material 

Poly  1  or  Poly  2 

Plate  Width 

50  jum 

Plate  Length 

100  /xm 

Support  Beam  Length 

60  /tm 

Support  Beam  Width 

2  /tm 

Number  of  Combs 

24 

Comb  Overlap 

33  /tm 

Comb  Length 

35  /tm 

Torsion  Bar  Width 

4  /tm,  3  /tm,  2  /tm 

Torsion  Bar  Length 

25  /tm,  10  /tm  (4  /tm  width  only) 

frequencies  were  calculated  using  a  MATLAB  script  file,  titled  ksys.m;  the  code  is  provided 
in  Appendix  C.  The  predicted  resonant  frequencies  are  listed  in  Table  4.6  with  each  of  the 
parameters  used  to  calculate  the  resonant  frequency.  In  Table  4.6,  ksys  is  the  system  spring 


Table  4.6  Comb  1  gyroscope  parameters. 


Parameter 

Value  by  structural  layer 

Poly  1 

Poly  2 

ksys 

2.192 

2.211 

M 

2.342  X  10"“ 

1.713  X  10-“ 

Msup 

5.620  X  10-13 

4.112  X  10-13 

Fr 

48.47  X  103 

56.92  X  103 

constant  in  Kgls"^,  M  is  the  mass  of  the  plate  and  any  combs  attached  in  grams,  Msup  is 
the  mass  of  the  support  springs  in  grams,  and  Fr  is  the  device  resonant  frequency  in  Hertz. 
In  all  of  the  comb  gyroscope  designs,  an  array  of  holes  may  be  seen.  These  holes  ensure 
a  complete  etch  of  oxide  from  underneath  the  large  plates,  and  they  reduce  the  viscous 
damping  in  the  output  oscillation  mode.  The  second  design  approach  to  the  tuning  fork 
gyroscope,  named  comb  2,  uses  a  large  plate  capacitor  as  the  sense  mechanism.  Figure  4.4 
shows  a  fabricated  device.  Variations  of  this  type  of  gyroscope  are  listed  in  Table  4.7.  As 


Figure  4.4  Comb2  gyroscope  design  In  MUMPS  7. 
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Table  4.7  Variations  and  Geometry  of  the  Comb  2  Gyroscope  Design. 


Parameter 

Value(s) 

Structural  Material 

Poly  1  or  Poly  2 

Plate  Width 

50  pm 

Plate  Length 

100  pm 

Support  Beam  Length 

60  pm 

Support  Beam  Width 

2  pm 

Number  of  Combs 

12 

Comb  Overlap 

33  pm 

Comb  Length 

35  pm 

Torsion  Bar  Width 

4  pm,  3  pm,  2  pm 

Torsion  Bar  Length 

25  pm,  10  pm  (4  pm  width  only) 

seen  from  this  table,  the  variations  are  the  same  as  for  the  first  approach  and  are  similarly 
motivated.  The  results  of  running  the  ksys  routine  for  this  design  and  its  variations  are 
given  in  Table  4.8.  The  resonant  frequencies  from  the  two  different  designs  are  identical. 


Table  4.8  Comb  2  gyroscope  parameters. 


Parameter 

Value  by  structural  layer 

Poly  1 

Poly  2 

ksys 

2.192 

2.211 

M 

2.342  X  10-“ 

1.713  X  10-“ 

Msup 

5.620  X  10-^3 

4.112  X  10-^^ 

Fr 

48.47  X  10^ 

56.92  X  10^ 

The  additional  mass  contributed  by  the  fingers  in  the  first  design  is  not  large  enough  to 
contribute  to  the  effective  mass  of  the  plate. 


J^.2.2  Planar  Gyroscope.  For  the  MUMPS  7  fabrication  run,  the  errors  in  the 
planar  gyroscope  layout  were  corrected  and  the  device  was  implemented  in  two  main  forms. 
The  first  is  shown  in  Figure  4.5,  here  the  device  is  close  to  square.  The  device  specifics 
are  given  in  Table  4.9.  The  second  approach  used  a  more  rectangular  shape  and  is  shown 
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Figure  4.5  Planar  gyroscope  design  1  in  MUMPS  7. 


Table  4.9  Variations  and  Geometry  of  the  First  Planar  Gyroscope  Design. 


Parameter 

Value(s) 

Structural  Material 

Poly  2 

Structure  Dimensions 

320  /im  X  290  /im 

Drive  Plate  Width 

80  fj,m 

Drive  Plate  Length 

80  /im 

Outer  Gimbal  Width 

2  pm,  4  pm,  6  pm,  8  pm 

Outer  Gimbal  Length 

10  pm 

Inner  Gimbal  Dimensions 

2  pm  X  5  pm  and  6  pmxlSOpm 

Number  of  Sense  Combs 

10 

Comb  Overlap 

35  pm 

Comb  Spacing 

2  pm 

Poly  2  Mass  Area 

23  pm  X  56  pm 

Gold  Mass  Area 

17  pm  X  50  pm 

Mass  location  (D  and  R  for  torque  equations) 

D=  75  pm  R=  30  pm 

in  figure  4.6.  The  device  specifics  and  variations  are  given  in  Table  4.10.  The  rectangular 
approach  allows  for  a  larger  sustained  oscillation  amplitude  but  limits  the  size  of  the  proof 
masses.  The  more  square  the  design,  the  larger  the  proof  masses  may  be,  but  the  smaller 
the  amplitude  of  the  sustained  oscillation.  The  output  torque  is  directly  proportional  to 
both  the  mass  and  the  amplitude  of  the  sustained  oscillation.  In  both  designs,  the  device 
has  a  sustained  oscillation  about  the  external  gimbal  excited  through  vertical  electrostatic 
actuation.  The  output  sense  mechanism  of  each  design  is  also  identical:  two  differential 
comb  electrodes.  There  were  two  elements  varied  in  both  the  designs.  The  first  was  the 
width  of  the  outer  flexure.  The  width  of  the  flexure  is  varied  to  adjust  the  device  resonant 
frequency.  The  second  used  two  different  proof  mass  placements  in  order  to  attempt  to 
maximize  the  output  torque,  the  first  placed  the  proof  masses  on  the  output  axis  and  the 
second  placed  the  proof  masses  off  the  output  axis.  Figures  4.7  and  4.8  show  the  two 
mass  placements.  The  first  mass  placement  would  result  in  zero  output  torque  about 
the  output  axis.  As  a  result,  the  devices  with  the  mass  placed  on  the  output  axis  will 
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Figure  4.7  Planar  gyroscope  mass  placement  on  the  output  axis  in  MUMPS  7. 
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Figure  4.8  Planar  gyroscope  mass  placement  olf  the  output  axis  in  MUMPS  7. 
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Table  4.10  Variations  and  Geometry  of  the  Second  Planar  Gyroscope  Design. 


Parameter 

Value(s) 

Structural  Material 

Poly  2 

Structure  Dimensions 

320  pm  X  210  ^m 

Drive  Plate  Width 

50  pm 

Drive  Plate  Length 

70  pm 

Outer  Gimbal  Width 

2  pm,  4  pm,  6  pm,  8  pm 

Outer  Gimbal  Length 

10  pm 

Inner  Gimbal  Dimensions 

2  pm  X  5  pm  and  6  /imx90/im 

Number  of  Sense  Combs 

10 

Comb  Overlap 

35  pm 

Comb  Spacing 

2  pm 

Poly  2  Mass  Area 

20  pm  X  40  pm 

Gold  Mass  Area 

14  pm  X  35  pm 

Mass  location  (D  and  R  for  torque  equations) 

D=  56  pm  R=  20  pm 

not  function.  The  second  placement  used  two  dumbell-like  structures.  Examination  of  the 
output  torque  equations  suggests  that  a  single  dumbell  structure  located  at  the  center  of 
the  output  axis  would  produce  a  larger  output  signal.  The  distance  out  from  the  center 
of  the  gyroscope  generates  a  negative  term  in  the  final  output  torque  equations.  If  the 
dumbell  structure  is  located  at  the  center  this  term  is  zero,  it  contributes  nothing  to  the 
output  torque. 


4-2.3  Ring  Gyroscope.  The  third  type  of  gyroscope  investigated  was  the  ring 
gyroscope.  This  gyroscope  was  first  investigated  at  the  University  of  Michigan  [10],  where 
the  device  was  constructed  from  electro-formed  nickel  yielding  a  high  aspect  ratio.  The 
aspect  ratios  attainable  in  the  MUMPS  process  make  the  Michigan  approach  to  driving 
and  sensing  the  ring’s  motion  impractical.  A  comb  drive  was  developed  to  maximize  the 
interface  area  when  driving  the  ring.  Two  methods  of  sensing  the  output  motion  were 
developed.  The  first  was  to  use  a  set  of  combs  identical  to  the  drive  combs  to  sense  the 
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output.  The  second  used  a  series  of  plates  suspended  over  poly  zero,  resulting  in  a  series 
of  capacitive  transducers  in  parallel.  Figure  4.9  shows  the  design  implemented  using  drive 
and  sense  combs.  The  second  type  of  ring  gyroscope  is  shown  in  Figure  4.10.  In  these 


Figure  4.9  Ring  gyroscope  design  using  sense  combs. 


designs,  a  single  drive  electrode  has  an  initial  capacitance  of  0.764  fF  and  a  drive  electrode 
from  Putty’s  gyroscope  has  a  stationary  capacitance  of  0.737  fF.  The  drive  capacitance 
of  the  two  structures  is  almost  identical;  as  a  result  the  force  generated  by  the  two  drive 


Figure  4.10  Ring  gyroscope  using  capacitive  plates  to  sense  output  motion. 


electrodes  is  also  very  close.  The  sense  electrodes  for  the  first  ring  design  are  identical 
to  the  drive  electrodes.  This  is  also  the  case  for  the  Putty  gyroscope.  When  plates  are 
used  instead  to  sense  the  output,  the  stationary  capacitance  is  0.850  fF,  an  increase  over 
the  stationary  sense  capacitance.  The  capacitance  values  for  the  sense  electrodes  provide  a 
means  of  comparison  between  the  two  electrode  types.  These  electrodes,  however,  are  used 
to  sense  the  change  in  capacitance.  For  a  larger  stationary  capacitance  in  both  designs 
there  will  be  a  larger  change  in  capacitance  for  an  output  motion. 

Each  of  the  ring  structures  fabricated  on  the  MUMPS  7  fabrication  run  has  an 
inner  radius  of  150  fim  supported  by  eight  semicircular  support  springs.  In  addition 
eight  electrodes  surround  the  structure.  Variations  of  this  design  are  changes  in  the  sense 
mechanism  and  support  structure.  Two  sense  mechanisms  were  developed.  The  first  is  the 
use  of  combs  identical  to  the  drive  combs  and  the  second  is  the  use  of  plates  suspended 
over  poly  zero  in  order  to  form  a  variable  capacitor.  The  nominal  comb  length  is  Z2fxm 
with  a  5fim  spacing  this  spacing  resulted  from  the  choice  of  ring  diameter.  There  are  seven 
plates,  each  having  an  initial  overlap  area  of  240iJ,m.  Three  different  base  designs  were 
investigated.  The  first  and  second  are  cylinders,  the  first  with  a  diameter  of  10/rrn  and  the 
second  with  a  diameter  of  20fj,m.  The  third  base  design  is  a  star  shape;  two  squares  with 
a  side  length  of  20/xm  are  anchored  to  poly  zero,  where  one  square  is  rotated  45  degrees 
relative  to  the  other.  The  goal  of  the  star  base  is  to  make  the  support  more  rigid,  shifting 
non-critical  resonant  frequencies  higher  than  the  resonant  frequency  of  the  two  flexural 
modes.  The  cylindrical  approach  was  first  used  in  Putty’s  working  designs.  Figures  4.11 
and  4.12  show  the  two  approaches  to  base  design.  The  ring  structures  fabricated  on  this 
run  contained  no  layout  errors.  The  probe  electrode  structure  used  to  make  connections  to 
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Figure  4.11  Ring  gyroscope  cylindrical  base  design. 
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Figure  4.12  Ring  gyroscope  star  base  design. 
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the  ring  gyroscopes  did,  however,  have  one  layout  error.  An  electrode  was  not  connected 
to  one  of  the  probe  pads.  This  was  corrected  in  the  MUMPS  8  fabrication  run. 


4.3  MUMPS  8 


The  statistics  from  this  run  are  given  in  Table  4.11.  The  largest  change  that  takes 


Table  4.11  MUMPS  8  Fabrication  Statistics. 


Film 

Thickness  (A) 

Stress  (MPa) 

Resistivity  (ohm-cm) 

poly  0 

4696 

8.0  (C) 

poly  1 

14.1  (C) 

poly  2 

10.3  (C) 

2.5E-3 

metal  (Cr/Au) 

166  (T) 

3.5E-6 

nitride 

5580 

133  (T) 

- 

oxide  1 

20362 

- 

- 

oxide  2 

6133 

- 

- 

place  on  this  fabrication  run  is  the  creation  of  stacked  structures.  A  stacked  structure 
increased  the  thickness  of  the  mechanical  layer.  This  results  in  higher  capacitance  values 
and  makes  it  possible  to  apply  a  larger  force  to  a  structure.  The  procedure  for  creating  a 
stacked  structure  follows.  The  layer  names  are  consistent  with  those  used  in  the  process 
description  found  in  Chapter  2,  Section  2.3.  To  create  a  stacked  structure,  first  poly  layer  is 
grown;  this  layer  is  then  patterned,  leaving  a  large  square  of  poly  1.  The  second  oxide  layer 
is  then  grown  and  the  plp2via  glass  cut  made  exposing  the  poly  1  plate.  Poly  2  is  grown 
and  patterned,  when  poly  2  is  patterned,  it  is  overetched  by  175  percent.  An  overetch  is 
the  duration  of  the  etch  past  the  minimum  required  to  cut  the  specific  layer.  Poly  2  is 
1.5  microns  thick  and  poly  one  is  2.0  microns  thick.  As  a  result,  the  overetch  of  poly  2  is 
enough  to  pattern  both  the  poly  1  and  poly  2  layers.  Figure  4.13  shows  each  of  the  steps  in 
creating  a  double  structure.  At  the  time  the  MUMPS  8  fabrication  run  was  sent  out,  the 
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Figure  4.13  Creation  of  double  structures  using  the  MUMPS  process:  (a)  Structure  after 
the  second  oxide  growth,  (b)  Structure  after  the  plp2via  glass  cut.  (c) 
Structure  after  the  poly  2  etch  is  completed,  (d)  Structure  after  it  is  released 
showing  the  anchored  poly  1  around  the  stacked  structure. 
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amount  of  the  poly  2  overetch  was  unknown;  as  a  result,  the  poly  1  plate  was  expanded 
beyond  the  plp2via  cut  and  anchored  to  prevent  excess  poly  from  floating  around  the  chip. 
Figure  4.14  shows  a  side  view  of  a  fabricated  double  structure.  Figure  4.14  also  shows  a 


Figure  4.14  Side  view  of  a  fabricated  double  structure. 


comparison  of  two  different  etch  mechanisms  used  at  MCNC.  The  first  etch  is  the  etch  that 
was  used  up  to  MUMPS  8,  shown  by  the  poly  1  in  the  right  of  Figure  4.14  (thinner  layer). 
The  sidewall  has  a  rough  appearance.  The  new  etch  mechanism  was  used  on  poly  2  and  is 
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shown  by  the  double  structure  (thicker  layer).  This  new  etch  yields  a  much  smoother  side 
wall  making  the  capacitance  equations  used  more  accurate. 


4.3.1  Vibrating  Comb  Gyroscope.  The  two  vibrating  comb  gyroscopes  developed 
for  the  MUMPS  7  fabrication  run  were  modified  for  the  MUMPS  8  run.  The  spring 
lengths  were  increased  to  reduce  the  resonant  frequency  of  the  devices  in  an  attempt  to 
increase  the  ability  to  drive  the  devices  at  their  resonant  frequencies.  The  comb  overlap 
was  increased  to  maximize  the  force  applied  by  the  drive  comb.  The  length  of  the  central 
support  spring  was  increased  to  maximize  output  sensitivity.  Comb  3  is  the  modified 
version  of  the  first  vibrating  comb  gyroscope  and  Comb  4  is  the  modified  version  of  the 
second  vibrating  gyroscope.  Table  4.12  has  the  design  parameters  for  the  Comb  3  design 
and  Table  4.13  contains  the  predicted  resonant  frequency  information.  Table  4.14  has  the 

Table  4.12  Variations  and  Geometry  Information  for  the  Comb  3  Gyroscope  Design. 


Parameter 

Value(s) 

Structural  Material 

Poly  1,  Poly  2  or  Stacked 

Plate  Width 

50  pm 

Plate  Length 

Support  Beam  Length 

80  pm 

Support  Beam  Width 

2  pm 

Number  of  Combs 

24 

Comb  Overlap 

30  pm 

Comb  Length 

40  pm 

Torsion  Bar  Width 

4  pm,  3  /rm,  2  ^m 

Torsion  Bar  Length 

40  pm 

design  parameters  for  the  Comb  4  design  and  Table  4.15  contains  the  predicted  resonant 
frequency  information. 
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Table  4.13  Comb  3  gyroscope  parameters. 


Parameter 

Value  by  structural  layer 

Poly  1 

Poly  2 

Stacked 

ksys 

2.809 

1.607 

5.003 

M 

2.305  X  10-^1 

1.800  X  10-11 

4.105  X  10-11 

Msup 

7.374  X  10-13 

5.760  X  10-13 

1.313  X  10-13 

Fr 

Table  4.14  Variations  and  Geometry  of  the  Comb  4  Gyroscope  Design. 


Parameter 

Value(s) 

Structural  Material 

Poly  1,  Poly  2,  or  Stacked 

Plate  Width 

50  /rm 

Plate  Length 

100  (xm 

Support  Beam  Length 

80  /im 

Support  Beam  Width 

2  nm 

Number  of  Combs 

12 

Comb  Overlap 

30  jum 

Comb  Length 

40  fim 

Torsion  Bar  Width 

4  //m,  3  /im,  2  /im 

Torsion  Bar  Length 

40  /im 

Table  4.15  Comb  4  gyroscope  parameters. 


Parameter 

Value  by  structural  layer 

Poly  1 

Poly  2 

Stacked 

ksys 

2.809 

1.607 

5.003 

M 

2.305  X  10-11 

1.800  X  10-11 

4.105  X  10-11 

Msup 

7.374  X  10-13 

5.760  X  10-13 

1.313  X  10-13 

Fr 

55.24  X  103 

47.27  X  103 

55.24  X  103 
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4.3.2  Planar  Gyroscope  and  Ring  Gyroscope.  The  planar  gyroscope  was  not 
changed  on  this  fabrication  run.  The  ring  gyroscope  design  was  not  modified.  A  third 
ring  gyroscope  was  created  in  order  to  take  advantage  of  the  double  structure  concept.  By 
doing  so  the  drive  electrode  capacitance  increases  to  1.16  fF.  This  increases  the  ability  to 
drive  the  structure. 


4.4  MUMPS  9 

The  MUMPS  9  fabrication  run  is  the  last  fabrication  run  used  in  this  thesis.  The 
statistical  data  on  the  run  is  given  in  Table  4.16.  There  were  changes  to  all  three  of  the 


Table  4.16  MUMPS  9  Fabrication  Statistics. 


Film 

Thickness  (A) 

Stress  (MPa) 

Resistivity  (ohm-cm) 

poly  0 

6243 

8.0  (C) 

1.5E-3 

poly  1 

19948 

3(C) 

2.2E-3 

poly  2 

15048 

6(C) 

3.0E-3 

metal  (Cr/Au) 

4923 

5(T) 

3.0E-6 

nitride 

6243 

67  (T) 

- 

oxide  1 

20008 

- 

- 

oxide  2 

7598 

- 

- 

gyroscopes  on  this  run  to  improve  on  the  previous  designs.  These  changes  are  outlined 
below. 


4.4-1  Vibrating  Comb  Gyroscope.  Testing  of  the  MUMPS  7  and  MUMPS  8  die 
illustrated  a  need  to  lower  the  spring  constant  for  the  vibrating  comb  gyroscope.  As  a 
result  the  support  springs  were  lengthened  in  both  the  comb  3  and  comb  4  designs.  The 
new  designs  are  named  comb  5  and  comb  6  respectively.  Tables  4.17  and  Table  4.18  list 
the  modifications  to  comb  3  and  comb  4  respectively. 
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Table  4.17  Variations  and  Geometry  Information  for  the  Comb  5  Gyroscope  Design. 


Parameter 

Value(s) 

Structural  Material 

Poly  2  or  Stacked 

Plate  Width 

50  jum 

Plate  Length 

100  /xm 

Support  Beam  Length 

200  fim 

Support  Beam  Width 

2  /tm 

Number  of  Combs 

24 

Comb  Overlap 

30  /xm 

Comb  Length 

40  /xm 

Torsion  Bar  Width 

4  /xm,  3  /xm,  2  /xm 

Torsion  Bar  Length 

40  /xm 

Table  4.18  Variations  and  Geometry  of  the  Comb  6  Gyroscope  Design. 


Parameter 

Value(s) 

Structural  Material 

Poly  2,  or  Stacked 

Plate  Width 

50  /xm 

Plate  Length 

100  /xm 

Support  Beam  Length 

200  /xm 

Support  Beam  Width 

2  /xm 

Number  of  Combs 

12 

Comb  Overlap 

30  /xm 

Comb  Length 

Torsion  Bar  Width 

4  /xm,  3  /xm,  2  /xm 

Torsion  Bar  Length 

40  /xm 
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The  largest  concern  with  these  designs  is  that  the  poly  may  sag  and  stick  to  the 


substrate. 


J^.J^.2  Planar  Gyroscope.  The  planar  gyroscope  design  was  changed  to  contain 
a  single  dumbell  like  structure  located  at  the  center  of  the  inner  axis.  This  increases  the 
output  torque  generated  by  a  given  input  rotation.  The  device  specifics  are  given  in  Table 
4.19.  Figure  4.15  shows  a  SEM  photograph  of  the  new  planar  gyroscope  design.  This  device 

Table  4.19  Variations  and  Geometry  of  the  First  Planar  Gyroscope  Design. 


Parameter 

Value(s) 

Structural  Material 

Poly  2 

Structure  Dimensions 

220  pm  X  320  pm 

Drive  Plate  Width 

30  pm 

Drive  Plate  Length 

85  pm 

Outer  Gimbal  Width 

2  pm 

Outer  Gimbal  Length 

10  pm 

Inner  Gimbal  Dimensions 

2  pm  X  5  pm  and  4  /imx290/im 

Poly  2  Mass  Area 

50  pm  X  220  pm 

Mass  location  (distance  from  center  axis) 

R=  50  pm 

will  maximize  the  output  torque  generated  for  a  given  input,  while  maintaining  symmetry. 
In  addition  the  interface  capacitance  of  this  type  of  planar  gyroscope  is  much  larger  than 
for  the  previous  designs.  The  larger  interface  capacitance  should  increase  device  sensitivity. 


4.4-3  Ring  Gyroscope.  There  were  two  main  modifications  to  the  ring  gyroscope 
on  this  fabrication  run.  First,  the  diameter  of  the  ring  was  doubled.  This  change  was 
dictated  when  initial  tests  failed  to  identify  the  device  resonant  frequency.  The  increase 
seems  dramatic;  in  the  initial  designs,  a  diameter  of  400  pm  was  used  in  order  to  prevent 
sagging  of  the  suspended  structure.  The  designs  developed  at  Michigan,  however,  used  a 
much  larger  diameter  of  1  mm  for  their  designs.  The  large  diameter  was  made  possible 
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Figure  4.15  Planar  gyroscope  4  design. 


largely  because  of  the  high  vertical  sidewalls.  The  new  ring  diameter  of  800  //m  in  poly 
may  result  in  the  structure  sagging  down  onto  the  substrate.  The  second  change  is  a  result 
of  an  ambiguity  created  in  the  first  designs.  In  the  ring  gyroscopes  on  the  MUMPS  7 
and  MUMPS  8  designs,  a  structure  with  eight  electrodes  was  using  and  each  electrode 
accounted  for  22.5  degrees  of  the  ring.  In  normal  device  operation,  a  resonant  node,  or 
maximum  of  the  resonant  pattern,  is  established  in  the  middle  of  the  sense  electrode.  For 
an  input  rotation,  that  node  moves  counter-clockwise  for  a  positive  rotation  and  clockwise 
for  a  negative  rotation.  In  both  the  devices  in  which  sense  combs  are  used  and  those 
in  which  sense  plates  are  used,  the  direction  of  rotation  will  be  unknown  until  the  node 
moves  11.25  degrees  in  either  direction.  This  severely  limits  the  device  sensitivity.  In  order 
to  reduce  the  effects  of  this  ambiguity,  the  sense  electrodes  have  been  made  smaller.  A 
single  sense  electrode  in  this  design  accounts  for  5  degrees  of  the  ring  and  the  resulting 
ambiguity  is  reduced  to  2.5  degrees.  The  cost  of  reducing  the  ambiguity  is  the  reduction 
in  the  nominal  capacitance  of  the  sense  electrodes.  The  ring  shown  in  Figure  4.16  is  the 
result  of  adjusting  for  these  concerns.  Electrodes  not  used  for  sensing  or  driving  the  ring 
will  be  grounded  with  respect  to  the  ring.  The  geometry  specifics  of  the  ring  gyroscope 
are  given  in  Table  4.20. 

Table  4.20  Variations  and  Geometry  of  the  Fourth  Comb  Drive  Gyroscope  Design. 


Parameter 

Value(s) 

Structural  Material 

Poly  2,  or  Stacked 

Diameter 

800  ^m 

Number  of  Support  Springs 

8 

Anchor  Type 

Star  Base 

Nominal  Comb  Overlap 

47  fim 

Nominal  Comb  Length 

56  fim 
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Figure  4.16  Modified  ring  gyroscope  design. 
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In  the  previous  four  sections,  design  considerations  for  the  three  types  of  gyroscopes 
were  discussed.  Specific  geometry  data  was  given  and,  where  applicable,  predicted  reso¬ 
nant  frequencies  were  generated.  The  next  section  covers  the  experimental  approach  and 
procedures  used  in  testing  the  devices,  including  the  procedures  used  to  process  the  die 
once  returned  by  MCNC. 

4.5  Experimental  Approach  /  Procedures 

In  this  section,  the  experimental  approach  and  procedures  are  covered.  First  a  general 
plan  of  attack  is  developed.  This  plan  is  followed  by  the  experimental  procedures  as 
implemented.  The  larger  task  of  testing  the  vibrating  gyroscopes  may  be  broken  down  into 
a  series  of  smaller  tests.  The  general  approach  to  testing  each  of  the  vibrating  gyroscope 
designs  is  given  below; 

1.  Attempt  to  drive  the  devices  to  resonance;  measure  the  corresponding  resonant  fre¬ 
quency  and  record  the  results. 

2.  Force  a  motion  identical  to  the  output  motion  physically  to  see  if  it  may  be  sensed. 

3.  Wire  bond  devices  out  and  verify  1  and  2. 

4.  Place  devices  into  a  rotating  field  to  apply  an  input  rotation  and  measure  output. 

In  the  following  sections  specific  procedures  used  will  be  provided.  An  equipment 
list  precedes  each  set  of  procedures. 

4.5.1  Release  Procedures.  When  the  chips  are  delivered  to  AFIT,  they  are  not 
yet  released.  The  release  process  involves  three  steps;  the  first  removes  any  photo-resist 
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on  the  die,  the  second  performs  the  silicon  dioxide  etch,  the  third  attempts  to  dry  the  die 
without  causing  the  devices  to  stick.  These  procedures  must  be  used  on  every  die  before 
it  may  be  tested. 

Equipment 

1.  1  Small  pail 

2.  4  150  ml  Beakers 

3.  1  Petri  dish 

4.  1  Hot  plate 

5.  1  Magnetic  stirrer  with  stir  bar 

6.  1  Tweezers 

7.  Rubber  gloves 

8.  Apron  and  Lab  Coat 

9.  Protective  Eye- ware 

Chemicals 

1.  Hydrofluoric  Acid  (HE)  49% 

2.  2-Propanol 

3.  Methanol 

4.  Acetone 

5.  De-ionized  (DI)  Water 
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Procedures 


1.  Clean  all  glassware  with  DI  water  and  blow  dry  with  nitrogen. 

2.  Fill  pail  about  half  full  of  DI  water.  Place  stir  bar  in  pail  off  to  one  side  so  die  may  be 
placed  in  pail  out  of  the  stir  bar’s  path.  Adjust  stirrer  speed  so  water  is  continuously 
moving. 

3.  Place  petri  dish  on  hot  plate.  Adjust  heat  so  the  dish  is  warm  to  the  touch. 

4.  Label  each  of  the  four  beakers  with  one  chemical  name  (HF,  Acetone,  2-Propanol, 
and  Methanol)  and  pour  approx.  50ml  of  the  appropriate  chemical  into  the  labeled 
beakers. 

5.  Place  the  die  to  be  released  into  the  acetone,  stirring  occasionally  for  5  minutes.  The 
acetone  is  intended  to  wash  any  excess  resist  away. 

6.  Remove  die  from  the  acetone  and  place  into  the  methanol,  stirring  occasionally  for 
3  minutes.  The  methanol  should  rinse  excess  acetone  and  resist  off  the  die. 

7.  Place  die  into  DI  water  bath  for  3  minutes  to  remove  the  methanol. 

8.  Place  die  into  HF  for  120  seconds.  The  HF  etches  the  silicon  dioxide,  releasing  the 
structure. 

9.  Remove  the  die  from  the  HF,  keeping  a  drop  of  liquid  on  top  of  the  die,  to  prevent 
stiction,  and  place  in  the  DI  water  bath  for  5  minutes. 

10.  Remove  the  die  from  the  DI  water  bath,  again  keeping  a  drop  on  top  of  the  die  and 
place  into  the  2-Propanol.  Leave  in  2-Propanol  for  5  minutes. 
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11.  Remove  die  from  2-Propanol  and  place  in  the  petri  dish.  If  the  plate  is  too  hot 
the  Propanol  will  boil,  flipping  the  chip.  This  should  be  avoided.  If  unsure  of  the 
temperature,  place  a  drop  of  Propanol  into  the  dish  and  watch  to  see  if  it  boils.  If  it 
does,  turn  the  heat  down  and  repeat.  The  goal  is  to  evaporate  the  Propanol  rapidly 
but  not  to  boil  it  off. 

12.  Inspect  the  die  under  a  microscope  to  verify  the  release.  If  the  die  is  not  completely 
released  return  the  die  to  the  HF  bath  and  repeat  steps  9-11. 

13.  If  the  die  has  been  released,  dispose  of  chemical  in  properly  marked  waste  bottles. 
The  DI  water  may  be  dumped  in  the  sink.  Rinse  and  dry  equipment  used. 

Steps  5-7  of  the  procedures  are  used  to  remove  any  excess  resist  on  the  die.  Steps  8  and 
9  perform  the  actual  release,  and  steps  10  and  11  attempt  to  dry  the  die  without  the 
devices  sticking  to  the  substrate.  When  a  device  sticks  to  the  substrate,  it  is  referred  to  as 
stiction.  Stiction  often  results  because,  as  the  device  dries,  the  meniscus  pulls  the  device 
to  the  substrate.  The  force  generated  is  a  function  of  the  surface  tension  of  the  liquid. 
Water  has  a  high  surface  tension  in  comparison  to  2-Propanol;  as  a  result,  if  the  water  is 
sublimated  with  2-Propanol,  there  is  a  better  chance  stiction  will  be  avoided.  If  the  devices 
do  stick  to  the  substrate,  it  is  not  possible  to  free  them.  Once  the  chips  are  released,  they 
may  be  tested. 

4.5.2  Resonant  Frequency  Testing.  The  resonant  frequency  testing  proved  to  be 
difficult.  Experimental  test  procedures  were  verified  using  a  comb  drive  lateral  resonator 
included  on  every  die  as  a  test  structure  by  MCNC  (see  Figure  4.17).  The  resonant 
frequency  of  this  test  structure  was  provided  by  MCNC  for  the  MUMPS  6  fabrication  run. 
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Three  methods  of  testing  for  resonance  were  investigated.  The  first  method  was  optical. 


Figure  4.17  MCNC  lateral  comb  resonator  test  cell. 


The  bandwith  of  resonance  is  very  small,  so  if  the  output  motion  is  visible,  resonance 
may  be  identified  by  slowly  varying  the  frequency  of  the  input  waveform  until  resonance 
is  observed.  The  second  method  used  a  network  analyzer  to  measure  resonance.  The 
third  and  final  method  involved  monitoring  the  current  output  by  the  voltage  source;  at 
resonance,  the  current  should  increase  as  more  energy  is  coupled  from  the  electrical  to  the 
mechanical  system.  Specific  procedures  for  testing  the  devices  are  included  in  this  section. 
An  equipment  list  precedes  each  of  the  procedure  sections. 


Each  of  the  gyroscopes  has  a  stationary  electrode  and  a  free,  or  moving  electrode. 
For  the  vibrating  comb  and  ring  gyroscopes,  these  are  a  set  of  combs.  On  the  planar 
gyroscopes,  the  stationary  electrode  is  a  poly  0  pad  and  the  moving  electrode  is  the  poly 
2  plate. 

4. 5. 2.1  Optical  Identification.  Previous  work  on  lateral  comb  resonators 
indicated  that  the  resonant  frequency  of  structures  may  be  determined  by  visual  identifi¬ 
cation  [27].  The  first  tests  were  conducted  based  on  this  approach.  The  procedures  used 
to  identify  the  resonant  frequency  visually  are  given  below. 

Equipment 

1.  High  voltage  DC  power  supply 

2.  HP  function  generator 

3.  Micromanipulator  Probe  Station  with  2  probes 

Procedures 

1.  Connect  the  stationary  electrode  to  the  output  of  the  function  generator  using  a 
BNC  to  BNC  cable.  Disable  the  function  generator  output  to  protect  devices  from 
damage,  adjust  the  amplitude  to  20  V  peak  to  peak,  and  set  the  initial  frequency  to 
10  kHz  (well  below  the  frequency  of  interest). 

2.  Connect  the  second  probe  from  the  DC  power  supply  to  the  structure.  Place  a 
DC  bias  on  the  structure  to  facilitate  resonance.  The  DC  bias  acts  as  a  static 
load,  increasing  the  magnitude  of  the  resonance  spike.  An  initial  value  of  60  V  was 
attempted.  This  value  was  increased  to  150  Vdc  for  testing  of  the  fabricated  devices. 
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3.  Enable  the  output  of  the  function  generator  and  slowly  begin  raising  the  frequency 
of  the  output  waveform.  At  resonance,  the  device  motion  will  be  visible. 

4.  Once  motion  is  visible,  it  is  possible  to  make  small  adjustments  to  the  function 
generator’s  output  frequency  so  the  maximum  displacement  is  observed. 

The  above  method  is  a  not  as  accurate  as  other  methods.  Because  the  bandwith  of 
resonance  is  small,  it  does  provide  a  good  first  measure  of  at  the  device  resonant  frequency. 

4. 5. 2. 2  Network  Analyzer  Measurements.  The  network  analyzer  provides 
a  reliable  method  of  measuring  the  resonant  frequency.  Several  experimental  set-ups  were 
used  to  conduct  resonant  frequency  measurements  on  the  vibrating  comb  and  vibrating 
ring  gyroscopes  using  the  network  analyzer.  Each  method  is  listed  below,  along  with 
variations  attempted. 

Equipment 

1.  1  HP  4195A  Network  /  Spectrum  Analyzer 

2.  2  HP  41800A  Active  probes 

3.  Micromanipulator  probe  station  with  3  probes 

4.  HP-IB  compatible  printer 

5.  High  voltage  DC  power  supply 

Procedures 

1.  The  network  measurements  require  applying  a  source  input  to  one  fixed  electrode  and 
measuring  the  output  from  a  second.  For  the  vibrating  comb  gyroscope,  the  input 
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is  connected  to  the  center  structure  and  the  output  is  measured  on  one  of  the  two 
other  fixed  electrodes.  For  the  ring  gyroscope,  the  input  is  placed  on  one  comb  drive 
and  the  output  is  placed  on  the  fixed  electrode  180  degrees  from  the  input.  Figure 
4.18  shows  a  general  diagram  of  the  network  analyzer,  showing  functional  blocks  and 
labeling  the  input  and  output  connections. 


Soft  keys  Measure  function  keys 


Probe  Power  S2  output 

Figure  4.18  Network  analyzer  block  diagram  labeling  output  and  input  connections. 

2.  All  measurements  were  conducted  using  active  probes  in  order  to  minimize  system 
noise.  The  probes  are  connected  to  the  R1  and  T1  Network  Analyzer  inputs.  The 
power  of  one  active  probe  is  connected  to  the  probe  power  connection  next  to  the 
SI  output  (source  1  output)  and  the  second  active  probe’s  power  is  connected  to 
the  probe  power  port  to  the  right  of  the  S2  output.  The  SI  output  is  connected  to 
the  Rl  input  (reference  1  input)  using  a  BNC  to  BNC  cable  and  a  dual  male  single 
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female  BNC  adapter.  The  female  connection  of  the  adapter  is  connected  to  the  probe 
attached  to  the  input  fixed  electrode. 

3.  The  Tl  input  (test  channel  input)  is  connected  to  the  probe  attached  to  the  output 
fixed  electrode. 

4.  The  DC  bias  from  the  network  analyzer  is  not  sufficient  for  these  tests.  The  high  volt¬ 
age  power  supply  is  used,  and  is  attached  to  the  free  electrode  (suspended  structure) 
and  an  initial  voltage  of  GOV  is  applied. 

5.  The  procedures  for  calibrating  the  network  analyzer  for  a  through  transmission  mea¬ 
surement  are  found  on  pages  4-19  and  4-20  of  the  network  analyzer  user  manual  [35]. 
In  order  to  ensure  the  accuracy  of  the  calibration,  both  the  input  and  output  probes 
were  placed  on  the  same  pad  to  short  the  input  to  the  output. 

6.  The  amplitude  of  the  output  oscillation  on  SI  was  adjusted  to  its  maximum  value  of 
1.26  V. 

7.  Different  frequency  ranges  were  investigated  to  increase  the  resolution.  The  network 
analyzer  breaks  the  frequency  range  into  a  maximum  of  401  pieces,  making  the 
minimum  step  size  equal  to  the  frequency  span  divided  by  401.  A  range  of  2  kHz  was 
used  to  ensure  a  high  resolution.  The  2  kHz  span  was  adjusted  to  cover  frequencies 
from  15  kHz  to  75  kHz  for  the  comb  and  ring  gyroscopes. 

8.  No  output  was  seen  for  the  vibrating  comb  or  ring  gyroscopes.  As  a  result,  the  DC 
bias  was  increased  to  150  V  and  the  output  was  amplified  to  +/-  10  V  peak-to-peak 
sinusoid.  The  network  analyzer  was  recalibrated  with  the  amplifier  in  line,  and  the 
frequency  range  from  15  kHz  to  75  kHz  was  again  investigated. 
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9.  The  second  modification  was  to  raise  the  substrate  potential  to  the  DC  bias  level 
in  an  attempt  to  reduce  the  lateral  friction  on  the  structures.  This  was  attempted 
both  with  and  without  the  amplifier  in  line.  The  amplifier  eventually  was  removed 
because  its  response  characteristics  dominated  the  measurements. 

10.  In  an  attempt  to  increase  the  quality  factor  of  the  resonance,  a  low  temperature 
probe  station  was  used  to  evacuate  the  chamber  to  ISmTorr.  The  measurements 
were  conducted  with  both  the  amplifier  in  and  out  of  line  with  the  source.  The 
measurements  were  conducted  in  a  vacuum  in  an  attempt  to  increase  the  magnitude 
of  the  resonant  spike. 

The  verification  and  results  of  these  experimental  procedures  are  provided  in  the 
following  chapter. 

4. 5. 2. 3  Current  Monitoring.  The  final  method  investigated  was  monitoring 
the  current  sourced  by  a  function  generator  as  the  input  frequency  was  swept.  At  mechan¬ 
ical  resonance,  the  current  sourced  should  increase  as  the  coupling  between  the  electrical 
and  mechanical  systems  increases.  The  maximum  current  observed  should  be  at  resonance. 

Equipment 

1.  Function  generator 

2.  Resistors  (1  ohm, 10  ohm, 100  ohm) 

3.  Volt  meter 

4.  Micromanipulator  probe  station  2  probes. 
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Procedures 


1.  Connect  the  function  generator  output  in  series  with  the  resistor  to  one  fixed  drive 
electrode. 

2.  Connect  the  volt  meter  across  the  resistor. 

3.  Apply  a  DC  bias  to  the  suspended  structure.  A  voltage  of  150  Vdc  is  used. 

4.  Adjust  the  output  of  the  function  generator  to  get  a  20  V  peak  to  peak  sinusoid  with 
an  initial  frequency  of  10  kHz. 

5.  Slowly  sweep  the  input  waveform  frequency  and  monitor  the  voltage  across  the  re¬ 
sistor. 

6.  Repeat  with  different  resistor  values  until  valid  output  is  obtained 

The  procedures  are  presented  in  the  order  that  they  were  attempted.  The  network 
measurements  taken  were  chosen  because  of  frequency  limitations  on  the  impedance  and 
reflectance  /  transmission  test  sets.  Reflectance  /  transmission  tests  sets,  used  to  determine 
the  two  port  network  parameters  (S-parameters),  have  a  minimum  frequency  threshold  of 
lOOkHz.  A  minimum  frequency  of  lOOkHz  is  approximately  two  times  the  frequencies  of 
interest.  The  impedance  test  sets  available  have  similar  frequency  limitations. 

4-6  Chapter  Summary 

This  chapter  presented  the  experimental  approach  and  procedures  used  in  completion 
of  this  thesis.  Device  descriptions  along  with  design  variations  were  presented  in  sections 
corresponding  to  the  fabrication  runs  the  designs  were  included  on.  The  procedures  used 
were  then  presented  along  with  equipment  lists  to  simplify  implementation.  The  following 
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chapter  will  provide  validation  of  the  procedures  and  the  observed  results  when  applied  to 
test  devices. 
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V.  Experimental  Validation  and  Results 


The  previous  chapter  outlined  the  general  approach  to  device  testing  and  then  de¬ 
scribed  the  specific  procedures  implemented.  This  chapter  provides  the  validation  of  those 
procedures  and  the  results  of  implementing  the  procedures  on  experimental  devices.  In 
addition,  a  discussion  of  the  observed  experimental  results  are  discussed. 

5. 1  Optical  Identification  of  Device  Resonant  Frequency 

Verification  of  the  optical  identification  procedures  was  conducted  on  die  from  the 
MUMPS  6  fabrication  run.  Predicted  resonant  frequencies  were  provided  by  MCNC  and 
six  different  MCNC  test  devices  were  investigated.  The  results  provided  from  MCNC 
were  based  on  a  50  volt  DC  bias  applied  to  the  suspended  structure  and  a  16  volt  peak 
to  peak  input  sinusoid  applied  to  the  drive  combs.  A  50  volt  DC  bias  was  also  used  to 
validate  the  experimental  procedures  outlined  in  Chapter  4;  however,  an  input  waveform 
of  an  18  volt  peak  to  peak  sinusoid  was  used.  The  amplitude  of  the  input  waveform  used 
by  MCNC  was  unknown  when  the  tests  were  conducted.  This  resulted  in  the  difference 
in  the  amplitude  of  the  waveform  used  to  conduct  the  tests  and  the  waveform  used  by 
MCNC.  The  results  provided  by  MCNC  indicate  an  output  oscillation  amplitude  of  3  pm  to 
5  /xm,  which  was  experimentally  verified.  The  results  of  testing  for  resonant  frequency  are 
provided  in  Table  5.1.  These  results  show  that  optical  identification  may  be  used  to  roughly 
determine  resonant  frequency  on  devices  for  which  the  output  motion  is  visible.  There  was 
no  output  oscillation  visible  below  25  kHz  or  above  28  kHz.  In  addition,  there  were  no 
significant  fabrication  errors  on  this  run.  A  larger  variation  in  resonant  frequencies  was 
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Table  5.1  Experimental  results  from  optical  resonant  frequency  identification  tests. 


Device  Name 

Observed  Resonant  Frequency 

MCNC 

26.3  kHz 

Device  1 

26.8  kHz 

Device  2 

26.9  kHz 

Device  3 

26.8  kHz 

Device  4 

26.6  kHz 

Device  5 

26.7  kHz 

Device  6 

26.6  kHz 

seen  on  the  MUMPS  7  and  MUMPS  8  fabrication  runs,  possibly  as  a  result  of  fabrication 
inconsistencies.  A  full  discussion  of  fabrication  problems  in  MUMPS  7  and  MUMPS  8  is 
included  at  the  end  of  this  chapter. 

The  optical  approach  to  resonant  frequency  determination  presents  a  simple  but 
crude  method  of  identifying  the  resonant  frequency.  Once  the  procedure  was  verified  it 
was  applied  to  the  vibrating  comb  and  ring  devices  fabricated  on  the  MUMPS  7,  MUMPS  8, 
and  MUMPS  9  runs.  Using  a  60  volt  DC  bias  and  a  36  volt  peak  to  peak  sinusoid,  the 
frequency  of  the  input  waveform  was  swept  slowly  from  15kHz  to  150kHz,  but  no  output 
motion  was  observed.  The  DC  bias  applied  to  the  structure  was  then  increased  to  100 
volts  and  then  followed  by  a  DC  bias  of  150  volts.  No  output  motion  was  observed. 


5.2  Network  Analyzer  Measurement  Validation  and  Resonant  Frequency  Results 

The  use  of  the  network  analyzer  was  validated  using  the  MCNC  lateral  comb  res¬ 
onators  on  the  MUMPS  7  die.  The  resonant  frequency,  of  the  MCNC  test  structures,  was 
first  measured  using  the  optical  method,  validated  in  the  previous  section.  The  network  an¬ 
alyzer  was  then  connected,  as  outlined  in  the  procedures  in  Chapter  4,  and  measurements 
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were  taken.  The  resonant  frequency  determined  optically  was  21.5  kHz  with  a  50  volt  DC 
bias  and  a  36  volt  peak-to-peak  sinusoid.  The  network  analyzer  measured  frequency  was 
21.3  kHz  also  at  a  50  volt  DC  bias.  The  output  of  the  network  analyzer  consists  of  two 
curves;  the  magnitude  of  the  test  channel  divided  by  the  magnitude  of  the  reference  channel 
(output  over  input)  in  mdB,  and  the  phase  in  degrees.  Several  different  DC  bias  voltages 
were  tested.  The  results  indicate  a  correlation  between  the  DC  bias  voltage  and  the  device 
resonant  frequency.  This  result  is  expected,  since  the  DC  bias  voltage  applies  a  static  load 
to  the  structure.  The  larger  this  load  is,  the  higher  the  resonant  frequency  will  be.  A  good 
analogy  to  this  effect  is  a  guitar  string.  The  resonant  frequency  of  the  string  is  adjusted 
by  changing  the  tension  on  the  string.  The  application  of  a  static  load  has  changed  the 
string’s  resonant  frequency.  Network  analyzer  outputs  are  shown  in  Figures  5. 1-5. 7  for  DC 
biases  from  0  volts  to  150  volts.  The  plot  for  the  DC  bias  of  50  volts  (Figure  5.3)  was  used 
to  compare  the  optical  results  and  the  network  analyzer  results  for  consistency.  In 
addition  to  increasing  the  resonant  frequency.  Figures  5.1-  5.7  show  that,  as  the  static  load 
is  increased,  the  resonant  spike  becomes  more  asymmetric.  The  dependence  of  resonant 
frequency  on  DC  bias  indicates  the  need  to  test  experimental  devices  at  a  set  of  different 
DC  biases.  The  initial  DC  bias  (60  volts)  was  increased  first  by  15  volts  DC  then  in  25  volt 
increments  to  its  maximum  value  of  150  volts.  Both  the  vibrating  comb  and  ring  gyroscope 
were  investigated.  Resonance  could  not  be  determined  for  experimental  vibrating  comb 
and  ring  devices  using  this  method.  Possible  explanations  for  the  inability  to  identify  the 
resonant  frequencies  of  experimental  devices  using  this  method  are  provided  in  Chapter  6. 

In  an  attempt  to  increase  the  amplitude  of  the  resonant  spike,  the  devices  were 
placed  in  a  vacnum,  chamber  at  15mTorr.  By  placing  the  devices  in  a  vacuum  the  effects 
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Figure  5.1  Network  analyzer  output  for  the  MCNC  comb  test  structure  with  a  zero  volt 
DC  bias  (MUMPS  7). 
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-58.70  44.00  T/R  -59.1170  dB 


DIV  DIV  START  16  300.000  Hz 

100.0m  2.000  STOP  26  300.000  Hz 

RBW:  30  Hz  ST:1.70  rain  RANGE:R“  10, T-  0dBm 

STEP-_25.001  HZ 


Figure  5.2  Network  analyzer  output  for  the  MCNC  comb  test  structure  with  a  25  volt 
DC  bias  (MUMPS  7). 
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NETWORK 

ft: REF  B:REF  o  mkR  22  125.233  Hz 

-58.60  44.00  T/R  -59.0532  dB 

[  dB  H  deg  1  6  20.1539  deg 


DIV  DIV  START  16  300.000  Hz 

100.0m  2.000  STOP  26  300.000  Hz 

RBW:  30  Hz  ST:1.70  min  RANGE:R-  10, T-  0dBm 

STEP-_25.001  HZ 


Figure  5.5  Network  analyzer  output  for  the  MCNC  comb  test  structure  with  a  100  volt 
DC  bias  (MUMPS  7). 


NETWORK 

ft --REF  B:REF  o  riKR  22  400.244  Hz 

-58.00  42.00  T/R  -58.937?  dB 


DIV  DIV  START  16  300.000  Hz 

200.0m  2.000  STOP  26  300.000  Hz 

RBW:  30  Hz  ST:1.70  min  RANGE:R-  10, T-  0dBm 

STEP-  25.001  HZ 


Figure  5.6  Network  analyzer  output  for  the  MCNC  comb  test  structure  with  a  125  volt 
DC  bias  (MUMPS  7). 
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DIV  DIV  START  16  300.000  Hz 

200.0m  5.000  STOP  26  300.000  Hz 

RBU:  30  Hz  ST:1.70  min  RANGE'.R-  10, T-  0dBm 

STEP-_25.001  HZ 

Figure  5.7  Network  analyzer  output  for  the  MCNC  comb  test  structure  with  a  150  volt 
DC  bias  (MUMPS  7). 

of  viscous  damping  decrease.  The  MCNC  test  structure  was  again  used  to  validate  the 
tests.  The  bias  voltage  required  for  an  output  response  to  be  observable  was  lowered 
significantly.  In  addition,  the  asymmetries  of  the  resonant  spike  increased  with  lower  DC 
bias  voltages.  The  results  are  shown  in  Figures  5.8-5.12.  Devices  on  both  MUMPS  7 
and  MUMPS  8  fabrication  runs  were  tested  in  the  vacuum  chamber.  Resonance  for  the 
experimental  devices  could  not  be  determined  using  this  method. 

5.3  Current  Monitoring  To  Identify  Resonance. 

The  current  monitoring  technique  to  identify  device  resonant  frequencies  has  not 
been  tested  in  this  thesis  and  is  left  as  a  possible  area  for  future  research. 
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Figure  5.8 


Figure  5.9 


g.pgp  ^  2^  070.000  Hz 

47.60  -30.00  T/R  46.4012  dB 


DIV  DIV  CENTER  21  000.000  Hz 

200.0m  2.000  SPAN  7  000.000  Hz 

RBW:  100  Hz  ST:40.2  sec  RANGE:R-  10, T-  10dBm 
0SC1-_1 .26E+00  V 


Network  analyzer  output  for  the  MCNC  comb  test  structure  with  a  0  volt  DC 
bias  in  a  vacuum  chamber  (MUMPS  7). 


NETWORK 

a: REP  BiREP  o  mkR 
46.60  -34.00  T/R 

1  dB  K  deg  1  6 


21  930.000  Hz 
46.268B  dB 

-39.0473  deg 


RBW:  30  Hz  ST =1.70  rain  RANGE :R- 

BIAS-_15.00  V 


10, T-  10dBm 


Network  analyzer  output  for  the  MCNC  comb  test  structure  with  a  15  volt 
DC  bias  in  a  vacuum  chamber  (MUMPS  7). 
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NETWORK 


A -.REF  B:REF  o  flKR  22  075.000  Hz 

46.80  -34.00  T/R  46.4010  dB 

[  dB  K  deg  1  B  -38.5229  deg 


DIV  DIV  START  20  000.000  Hz 

100.0m  1.000  STOP  24  000.000  Hz 

RBW:  30  Hz  ST:40.9  sec  RANGE:R-  10, T-  10dBm 

BIAS-  20.00  V 


Figure  5.10  Network  analyzer  output  for  the  MCNC  comb  test  structure  with  a  20  volt 
DC  bias  in  a  vacuum  chamber  (MUMPS  7). 


B:REF  o  MKR  23  075.000  Hz 

46.50  -34.00  T/R  46.0491  dB 

C  dB  ][  deg  1  e  -36.4525  deg 


DIV  DIV  START  20  000.000  Hz 

100.0m  1.000  STOP  24  000.000  Hz 

RBW.  30  Hz  ST:40.9  sec  RANGE:R-  10, T-  10dBm 
BIAS-  26.00  V 


Figure  5.11  Network  analyzer  output  for  the  MCNC  comb  test  structure  with  a  26  volt 
DC  bias  in  a  vacuum  chamber  (MUMPS  7). 
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24  390.000  Hz 

46.50  -34.00  T/R  45.85B1  dB 

[  dB  K  deg  3  8  -35.2635  deg 


DIV  DIV  START  31  000.000  Hz 

100. 0tn  1.000  STOP  25  000.000  Hz 

RBW:  30  Hz  ST;1.70  min  RANGE:R-  10. T-  10dBm 

BIAS-_30.00  V 


Figure  5.12  Network  analyzer  output  for  the  MCNC  comb  test  structure  with  a  30  volt 
DC  bias  in  a  vacuum  chamber  (MUMPS  7). 

5.4  Discussion  of  MUMPS  Fabrication  Process. 

Four  MUMPS  fabrication  runs  were  participated  in;  MUMPS  6,  MUMPS  7,  MUMPS  8, 
and  MUMPS  9.  Of  these  runs,  two  were  returned  without  fabrication  difficulties;  they  were 
the  MUMPS  6  and  MUMPS  9  fabrication  runs.  The  residual  stresses  of  the  material  layers 
in  these  runs  were  significantly  smaller  than  those  seen  in  the  MUMPS  7  and  MUMPS  8 
runs.  There  were  four  problems  with  the  MUMPS  7  and  MUMPS  8  fabrication  runs.  The 
first  problem  was  located  in  places  the  poly  1  layer  and  poly  2  layers  had  been  fused.  This 
was  not  a  significant  issue  for  the  devices  designed  as  resonators.  The  second  problem, 
however,  caused  significant  difficulty.  Where  the  topology  changes  the  polysilicon  layers 
will  be  thicker.  As  a  result,  if  the  polysilicon  etch  is  not  long  enough,  thin  lines  of  polysili¬ 
con,  called  stringers,  are  left  behind.  Figure  5.13  shows  how  stringers  are  created.  Once  the 
die  is  released  these  lines  stick  to  the  substrate.  For  electrostatic  devices,  these  thin  lines 
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Poly  1 
Oxide  1 


Reactive  Ion  Etch 


(a) 


(b) 

Figure  5.13  Stringer  creation  in  the  MUMPS  process:  (a)  Polysilicon  is  deposited  across 
the  wafer  and  then  etched  using  a  reactive  ion  etch.  The  polysilicon  layer 
is  thicker  where  contour  changes  occur,  requiring  a  longer  etch  time,  (b) 
Stringers  are  the  result  of  an  incomplete  reactive  ion  etch,  as  shown  in  this 
figure. 
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Figure  5.14  Stringers  shown  shorting  the  MCNC  test  structure  and  making  it  unusable. 

current  surge  is  sufficient  to  melt  the  device,  probe  tips,  and  often  create  a  significant 
depression  in  the  substrate.  While  stringers  were  present  on  the  MUMPS  7  fabrication 
run,  they  were  not  abundant  and  therefore  had  little  effect.  There  were  significantly  more 
stringers  on  MUMPS  8,  making  device  testing  difficult  and  the  results  questionable.  In 


of  polysilicon  or  stringers  align  themselves  with  the  electrostatic  field.  Often  the  result  is 
shorting  from  a  high  potential  (150V  DC)  to  a  low  potential  (OV  DC).  Figure  5.14  shows 
an  MCNC  test  structure  that,  because  of  stringers,  would  short  instantly.  The  resulting 
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addition  to  stringers,  the  residual  stress  of  the  mechanical  layers  increased  consistently 
from  MUMPS  6  through  MUMPS  8  and  then  decreased  again  in  MUMPS  9.  For  exam¬ 
ple,  the  residual  stress  of  the  poly  1  layer  increased  from  a  moderate  value  of  5.6  MPa 
on  MUMPS  6  to  approximately  2.5  times  that,  14.1  MPa,  on  MUMPS  8,  increasing  the 
device’s  resonant  frequency.  The  residual  stress  of  the  poly  1  layer  decreased  in  MUMPS  9 
to  3  MPa.  Variations  in  the  residual  stress  characteristics  led  to  large  variations  in  device 
resonant  frequencies.  In  addition,  the  lack  of  consistency  in  the  residual  stress  charac¬ 
teristics  makes  it  difficult  to  design  devices  with  specific  resonant  frequencies.  The  third 
problem  encountered  may  or  may  not  be  a  fabrication  error.  It  has  been  observed  that  a 
thin  film  is  present  over  many  poly  0  surfaces.  This  film  causes  devices  to  stick  in  posi¬ 
tion.  They  may  be  moved  laterally  with  probes;  however,  when  probe  pressure  is  removed, 
the  structure  stays  in  its  last  position.  The  film  causes  the  devices  to  stick  as  a  result 
of  increasing  the  friction  present.  Several  methods  were  investigated  to  remove  this  film, 
including  placing  the  die  in  both  organic  and  inorganic  solvents.  Neither  of  these  methods 
was  effective.  It  is  possible  that  the  film  was  the  result  of  contaminated  DI  water  or  a 
contamination  at  the  foundry.  While  the  majority  of  the  MCNC  test  structures  functioned 
on  the  MUMPS  7  die,  only  one  of  15  structures  tested  on  MUMPS  8  functioned  and  many 
others  were  not  tested  because  of  the  number  of  stringers  on  the  structures.  The  fourth 
challenge  presented  by  the  MCNC  fabrication  process  was  the  lead  time  from  sending  out 
devices  to  having  them  returned.  The  devices  sent  out  took  three  months  to  return  on 
each  of  the  MUMPS  6,  MUMPS  7,  and  MUMPS  8  fabrication  runs.  The  runs  were  also 
spaced  three  months  apart.  As  a  result,  fabricated  devices  were  returned  no  more  than 
two  or  three  days  before  the  next  run  was  due  to  be  sent  out.  There  was  no  time  to  test 
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devices  between  runs.  The  lead  time  has  been  decreased  and  was  two  and  a  half  months 


for  the  MUMPS  9  fabrication  run. 

5.5  Discussion  of  Resonant  Frequency  Testing  Results 

The  optical  and  network  analyzer  methods  of  testing  for  device  resonant  frequency 
was  conducted  on  devices  from  the  MUMPS  7,  MUMPS  8,  and  MUMPS  9  fabrication  runs. 
The  procedures  were  verified  using  the  MCNC  test  structure  and  have  proven  effective  in 
determining  the  resonant  characteristics  of  micro  mirrors  and  arrays  of  those  mirrors  [36] . 
There  are  two  possible  reasons  why  the  experimental  gyroscope  devices  failed  to  work. 
The  first  possibility  is  that  the  flexures  provided  too  large  a  spring  constant.  As  a  result, 
too  large  a  force  was  necessary  to  move  the  devices  both  at  atmospheric  pressure  and 
when  in  a  vacuum.  The  second  possibility  is  that  the  suspended  mass  was  large  enough 
to  cause  the  devices  to  sag  and  either  stick  to  the  substrate  or  to  move  the  free  drive 
combs  out  of  alignment  with  the  fixed  drive  combs.  Figure  5.15  shows  a  stacked  structure 
with  12  fingers.  The  device  plate  is  suspended  off  the  substrate  and  the  combs  are  still 
interdigitated.  Figure  5.16  shows  a  double  structure  with  24  combs,  all  other  geometry 
features  are  the  same.  The  picture  shows  the  combs  are  not  interdigitated  and  the  device 
plate  appears  to  be  resting  on  the  substrate.  Ten  to  fifteen  devices  of  each  type  were 
investigated,  all  showed  similar  results.  The  consistency  in  the  comb  alignment  indicates 
the  MUMPS  8  devices  have  one  of  three  problems.  The  first  possibility  is  the  four  flexures 
in  the  vibrating  members  are  too  long.  This  is  possible  as  the  spring  constant  in  the 
vertical  direction  is  identical  to  the  spring  constant  in  the  lateral  direction.  The  second 
possibility  is  the  length  of  the  output  flexure  may  be  too  large,  resulting  in  enough  sag 


5-14 


Figure  5.15  Stacked  comb  6  structure  showing  the  overlap  of  the  combs  (MUMPS  8). 
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to  lower  the  combs.  The  amount  of  drive  force  generated  is  decreased  linearly  with  the 
decreasing  overlap  area  of  the  combs.  The  third  possibility  is  that  stiction  caused  the 
structures  to  stick  to  the  substrate.  The  devices  were  located  in  different  quadrants  of 
the  die.  Stiction  may  have  been  a  larger  factor  in  one  quadrant  than  in  the  other.  The 
Comb  gyroscope  designs  developed  for  the  MUMPS  9  fabrication  run  were  based  on  a  high 
residual  stress  value.  The  MUMPS  9  run  ended  up  having  a  low  residual  stress  value.  As 
a  result,  the  devices  sagged  and  stuck  to  the  substrate. 

The  ring  gyroscopes  in  the  MUMPS  7  and  MUMPS  8  fabrication  runs  did  not  illus¬ 
trate  as  consistent  a  result.  Some  rings  were  suspended  off  the  substrate  as  expected  and 
others  were  not.  The  inconsistencies  indicate  the  cause  as  stiction.  Methods  have  been 
developed  to  minimize  the  effects  of  stiction  [30].  One  such  method  places  the  die  in  a 
pressure  chamber.  Any  liquid  on  the  die  is  sublimated  using  liquid  carbon  dioxide.  The 
chamber  is  then  brought  to  a  pressure  consistent  with  the  triple  point  of  carbon  dioxide. 
The  triple  point  refers  to  the  temperature  and  pressure  where  a  substance  may  exist  in 
the  liquid,  gas,  and  solid  phases.  The  carbon  dioxide  then  transitions  from  a  liquid  to 
a  gas  instantaneously.  Using  this  method,  a  100  percent  successful  release  rate  has  been 
reported  [30].  The  ring  gyroscope  design  included  on  MUMPS  9  had  twice  the  diameter. 
Again  the  choice  to  increase  the  diameter  was  based  on  the  high  residual  stress  of  the 
previous  two  runs.  The  lower  stress  in  MUMPS  9  caused  these  devices  to  sag  and  stick  to 
the  substrate. 

The  planar  gyroscope  designs  on  the  MUMPS  7  and  MUMPS  8  fabrication  runs  were 
tested  using  the  both  the  optical  and  network  analyzer  methods.  Resonance  could  not  be 
determined  applying  these  methods.  The  outer  flexures  on  these  designs  may  have  been 
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too  short  to  twist  properly.  The  MUMPS  9  planar  gyroscope  designs  were  based  on  a  high 
residual  stress.  With  the  lower  residual  stress  encountered  these  devices  sagged  and  stuck 
to  the  substrate. 

5. 6  Chapter  Summary 

This  chapter  presented  the  experimental  verification  of  two  methods  to  determine 
device  resonant  frequencies.  The  first  was  the  optical  method.  This  approach  leads  to  a 
good  initial  estimate  of  the  device  resonant  frequency.  The  second  used  a  network  analyzer 
to  measure  the  output  over  input  transfer  function.  This  method  provides  a  more  accurate 
measure  of  resonant  frequency.  Using  the  network  analyzer  a  dependence  of  resonant 
frequency  on  bias  voltage  was  observed  along  with  the  effects  of  testing  devices  in  a  vacuum. 
Both  the  optical  method  and  the  network  analyzer  methods  were  applied  to  experimental 
gyroscope  designs.  The  resonant  frequency  of  the  experimental  gyroscopes  could  not  be 
identified  using  either  method.  Possible  reasons  for  the  difficulty  in  determining  the  device 
resonant  frequency  and  suggested  solutions  are  provided  in  Chapter  6. 
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VI.  Conclusions  and  Recommendations 

This  chapter  presents  conclusions  and  recommendations  based  on  the  information  in 
previous  chapters.  The  conclusions  drawn  from  this  work  are  presented  first  followed  by 
recommendations  for  future  research  in  this  area. 

6. 1  Conclusions 

This  thesis  investigated  the  use  of  a  commercial  fabrication  process  to  fabricate  a 
micromachined  gyroscope  capable  of  being  integrated  into  an  INS  /  GPS  system  providing 
accurate  position  and  velocity  data.  The  advantages  of  the  micromachined  gyroscopes 
are  the  low  cost  and  small  size  of  the  devices.  Three  different  designs  were  investigated: 
the  vibrating  comb  gyroscope,  the  vibrating  ring  gyroscope,  and  the  planar  gyroscope. 
Operational  theory  was  developed,  from  first  principals,  for  the  vibrating  comb  and  pla¬ 
nar  gyroscopes.  The  ring  gyroscope  theory  was  initially  developed  at  the  University  of 
Michigan  [7].  Modifications,  however,  had  to  be  made  for  changes  in  the  drive  and  sense 
mechanisms. 

The  commercial  fabrication  process  available  was  MCNC’s  MUMPS  process.  Designs 
were  submitted  for  four  fabrication  runs  MUMPS  6-MUMPS  9.  Of  the  returned  test  die, 
the  MUMPS  6  and  MUMPS  9  runs  had  the  least  number  of  fabrication  errors  and  the  lowest 
residual  stresses.  There  was  considerable  variation  in  the  parameters  of  each  fabrication 
run,  making  it  difficult  to  design  devices  for  a  specific  resonant  frequency.  Stringers  were 
found  on  both  the  MUMPS  7  and  MUMPS  8  fabrication  runs.  The  stringers  presented  a 
larger  problem  on  the  MUMPS  8  run,  causing  electrical  shorts  that  resulted  in  destruction 
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of  devices  and  probes.  Before  a  commercial  fabrication  process  may  be  used  to  create  test 
devices  or  mass  produce  devices,  it  must  be  able  to  produce  consistent  reliable  results.  The 
MCNC  process  has  not  been  able  to  meet  these  criteria. 

Three  methods  were  used  in  an  attempt  to  evaluate  device  resonant  frequency:  op¬ 
tical  detection,  network  analyzer  measurements,  and  current  monitoring.  The  current 
monitoring  tests  have  not  been  completed  and  are  left  as  an  area  for  future  research.  The 
procedures  used  in  optical  detection  and  network  analyzer  measurements  were  verified 
through  the  use  of  a  test  comb  drive  resonator  structure  included  on  each  die  by  MCNC. 
The  procedures  were  then  applied  to  the  experimental  gyroscope  structures.  The  resonant 
frequencies  for  the  vibrating  comb  gyroscope,  the  vibrating  ring  gyroscope,  and  the  planar 
gyroscope  could  not  be  determined  using  either  of  the  first  two  methods.  There  are  several 
possible  explanations  for  the  inability  to  identify  the  resonant  frequency.  For  the  vibrating 
comb  gyroscopes: 

•  The  flexures  may  be  too  stiff;  as  a  result  too  much  energy  is  required  to  drive  the 
device. 

•  The  overlap  of  the  fixed  and  free  combs  may  be  reduced  enough  by  sag  to  increase 
the  energy  required  beyond  the  capabilities  of  the  equipment  available. 

•  Stiction  may  be  causing  devices  to  stick  to  the  substrate,  making  them  unusable. 

On  the  vibrating  ring  gyroscope,  the  flexures  may  be  too  thin,  causing  the  device  to  sag, 
resulting  in  stiction  or  in  moving  the  drive  combs  out  of  position,  reducing  the  applied 
force.  In  addition,  the  drive  force  is  applied  over  a  larger  arc  length  than  on  the  devices 
developed  at  the  University  of  Michigan.  As  a  result,  the  effect  of  applying  an  equal  force 
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is  less  on  devices  designed  for  this  thesis.  The  following  section  provides  recommendations 
for  continuing  research  in  this  area. 

6.2  Recommendations 

This  section  provides  a  set  of  recommendations  that  may  prove  helpful  in  continuing 
research  in  this  area. 

1.  Verify  and  attempt  resonant  frequency  identification  using  the  current  monitoring 
method  outlined  in  Chapter  4. 

2.  Investigate  other  possible  fabrication  sites. 

3.  Develop  test  structures  investigating  the  maximum  mass  that  can  be  suspended 
using  a  suspension  identical  to  that  of  the  comb  gyroscope.  The  devices  placed 
on  the  MUMPS  7  test  die  appear  to  be  too  stiff  and  devices  on  the  MUMPS  8  and 
MUMPS  9  test  die  appear  to  be  too  loose,  resulting  in  devices  sagging. 

4.  Develop  test  structures  investigating  the  effects  of  varying  the  support  springs  in 
the  ring  gyroscope.  Attempt  both  thinner  and  thicker  support  flexures;  observe  the 
effects  these  changes  have  on  the  comb  alignment. 

5.  Shorten  the  torsion  bar  length  of  the  comb  3  -  comb  6  gyroscopes.  The  device  sag 
may  be  a  result  of  too  much  mass  suspended  from  these  flexures. 

6.  The  overetch  of  poly  2  is  sufficient  to  make  a  poly  1,  poly  2,  poly  0  stacked  structure. 
Create  fixed  combs  from  there  three  layers  to  maximize  comb  interface  area  and,  as 
a  result,  maximize  the  applied  force. 
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7.  Alternative  flexure  designs  should  be  investigated.  Wander  springs  have  been  used 
with  success  in  lateral  resonators  [37].  The  effects  of  the  output  torque  on  this  type 
of  flexure  should  first  be  investigated. 

8.  Improvements  may  also  be  made  to  the  comb  drive  mechanism.  The  comb  drives 
as  implemented  result  in  a  vertical  as  well  as  horizontal  applied  force.  Methods 
have  been  developed  to  suppress  the  levitation  effect  and  attain  a  pure  horizontal 
force  [30]. 

9.  Attain  access  to  high  voltage  dc  supplies  and  ac  oscillators  to  maximize  the  energy 
applied  to  the  gyroscopes. 

10.  Initially  it  was  thought  the  devices  were  suspended  a  sufficient  height  off  the  sub¬ 
strate.  Tests  have  indicated  this  is  not  true  and  methods  of  preventing  stiction  must 
be  implemented.  Two  methods  have  been  used  successfully  on  previous  designs: 

(a)  Place  small  depressions  in  the  poly  2  structures  using  the  dimple  layer  to  mini¬ 
mize  the  area  that  results  in  stiction. 

(b)  Place  a  large  plane  under  the  suspended  structure  at  equal  potential. 

This  thesis  establishes  a  base  for  research  in  the  micromachined  gyroscope  area  using 
a  commercial  fabrication  process.  An  overview  of  three  device  designs  has  been  provided; 
in  continuing  this  research,  a  narrower  focus  is  necessary.  Many  opportunities  exist  to 
continue  this  research  at  all  levels  to  include:  Attempt  to  increase  the  mass  of  the  planar 
gyroscope  designs  to  increase  output  sensistivity,  re-design  of  the  flexure  systems  for  the 
comb  and  ring  gyroscopes,  and  alterations  to  the  comb  drive  mechanisms.  In  addition 
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effort  must  be  dedicated  to  developing  low  noise  sense  circuitry  and  a  means  of  interfacing 
that  circuitry  with  the  test  micromachined  structure. 
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Appendix  A.  Vibrating  Comb  Gyroscope  Output  Derivations 


In  this  appendix,  generic  equations  will  be  derived  for  the  compound  accelerations 
and  equations  for  output  torque  generated  for  Figure  A.l  In  Figure  A.l,  CML  and  CMR 


Figure  A.l  Block  diagram  for  generic  torque  derivations 

represent  the  left  and  right  centers  of  mass.  Position  vectors  to  CML  and  CMR  are  defined 
in  Equations  (A.l)  and  (A. 2).  The  input  rotation,  u;i„,  will  be  defined  by  Equation  (A. 3): 


Rl  =  -{jL  +  0.sin{u}t)xg  +  b^yg  -  H^Zg 


Rr  =  —{Vr  +  Qsin{u>t)Xg  +  brPg  —  HrZg 


= 
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(A.l) 

(A.2) 


(A.3) 


The  acceleration  that  results  from  an  input  rotation  is  calculated  using  the  same 
method  as  in  the  Coriolis  acceleration  example.  First  the  acceleration  experienced  by 
the  left  vibrating  member  is  Ccdculated.  The  calculation  for  the  acceleration  of  the  right 
vibrating  member  follows. 

Vl  =  [— Osm(wt)  —  QhJCOs{ijjt)  —  Hiujy  — 

+[-Hi(v^  +  [ri  +  9.sin{Ljt)]uj,]yg 

+[bL^x  +  [ri  +  ^sin{ut)]LOy]zg  (A.5) 

Ai  —  [[— O  +  [up'  +  Wy  —  ijol\'[l\sin{u>t)  —  [2f2w  +  Quj]cos{u>t) 

+[u^0Jy  -  oj,]bL  +  [Wj,  -  oj^oj,]Hl  +  [u^  -  oj^Prilxg 

+[[2c<;jli  +  +  uP}^]sin{ut)  +  2fluJ^wcos(wt) 

+  +  ^^1]hL  +  [<^y^z  —  +  [WxWy  +  ^z\l'L]yg 

+[[2u;yO  +  \u)y  +  (jj^Ux\‘^]sin{ix!t)  +  2Q,u}yU!cos{u!t) 

+  [(^X  +  +  ‘^y]rL]Zg  (A.7) 


The  acceleration  of  the  right  member’s  center  of  mass  is: 


Vr  =  +  Q.ujcos{(jjt)  —  H^LOy  —  brU^z]xg 

-  [Vr  +  9.sin{ij:t)]uj^]yg 

->r[Ki>^^  -  [r-r  +  0,sm{ijjt)]LVy\Zg 
- 

Ar  =  [[O  +  [— —  (jy  A  u)l]0]sin(u)t)  +  [2Clui  +  Qu}]cos(Lot) 

+  [u!^u>y  -  oj^]br  +  [-OJy  +  +  [-<^1  +  i^lVAXg 

+[— [2a;^J7  +  {ux<^y  +  —  2Q,uJzU>cos{Lot) 

+[^x  +  ‘^1%  +  Ko;^  -  iA:x\Hr  +  [-U>^UJy  -  !j,]rr]% 
A[—[‘iUytl  A  [d>y  A  ‘j^z'^x\'^]sin{u)t)  —  2Q.UyUicos{ut) 

A[CJx  +  Wj,Wj]6r  +  [-w^  +  WyliTr  +  [-OJxl^z  -  ojy]rr]Zg 

The  acceleration  of  the  mass  centers  causes  a  net  torque  about  the  output 
definition  of  torque  is  repeated  here: 


f  =  Rx  F 
F  =  mA 

The  torque  from  the  left  vibrating  member  is: 


T[ 


R\  X  Ff 


A-3 


(A.9) 

(A.IO) 


(A.ll) 

is.  The 

(A.12) 

(A.13) 


(A.14) 


T[  —  [bimL[[2u)y{l  +  [u}y  +  Lj^ujj;]Q]sin(ut)  +  2^UyUcos(ut) 

+  [Wj:  +  ^yl^z]bL  +  [“^x  +  +  [<^x^z  + 

+HLmL[[2u}2(l  +  (wx^y  +  (^z)H]sin{ut)  +  2flw^wcc»s(a;i) 
+  [ul  +  iJ^]bL  +  [l^yOJz  ~  +  ^^z]rL]]Xg 

+[(ri  +  ^lsin(LJt))mL[—[2(jjyQ  +  [d>y  +  —  2Q,UyUJCos{ujt) 

-[uj^  +  u}yLJ^]bL  -  -  [(z^xi^z  +  Wy]r-i] 

+HimL[[~Q  +  [u^  +  Wy  —  -  [2ilu>  +  ncj]cos(a;i) 

+  [ii>^Uy  -  a3j6i  +  [Wj,  -  -  u}'^,]rL]yg 

+[(ri  +  Q,sin(u!t))mL[—[2u)^Cl  +  {uxOJy  +  -  2Q.Lo^u>cos{u>t) 

-[W®  +  ^^z]^L  -  Kw^  -  UxIHl  -  [WxWy  +  w^]ri] 
-6i7ni[[- ii  +  [w^  +  ^^y  —  —  [2fiw  +  f2w]cos(wi) 

+  [WxWy  -  ^z]h  +  [Wy  “  +  [w^  “ 

The  torque  from  the  right  vibrating  member  is: 


n  =  mx  F? 

Tr  =  [brmr[—[2u)yQ,  F  [chj,  +  uj;s(jjj;]H]sin{ut)  —  2^ujyUcos{ut) 
+  [(hx  +  WyWj&r  +  [-ojI  +  cjl]Hr  +  [-OJ^OJ^  -  ujy]rr\ 
+  ifr^r[“[2w^fi  +  {^x‘zz>y  +  OJ z)^\sin{ujt)  —  2^iz)zUJC0s{u}t) 

+  [<^l  +  ^l]br  +  -  Wx]ffr  +  [-WxWj,  -  iz),]r^]]Xg 

+[[rr  +  ilsin{ut)]mr[-[2u}ytl  +  [uy  +  a;^ti;x]0]sm(a;t)  -  29.u)yU>cos{ut) 


(A.15) 


(A.16) 
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+  [u}^  +  CjJyU^]br  +  [-ul  +  u)l]Er  + 

+Hr'mr[[^  +  [-uj^  “  ‘*^y  +  >j^l]'^]sin{ut)  +  [itluj  +  ^u]cos{u)t) 

+  [u^UJy  -  uj,]br  +  [-W;,  +  (^a!U,]Hr  +  [-U}y  +  l^tVAyg 

+[[rr  +  fism(ut)]mr[— [2tJ^$7  +  —  2Cloj^u}cos{ut) 

+  [ul  +  +  [^y>^z  -  +  [-l^x^^y  “  ^^zVr] 

—briTirll^  +  [—oj^  ~  ^y  +  u)l]Q,]sin(uit)  +  [2$2u;  +  ^u]cos(Ljt) 

+  [uJxOJy  -  uj,]br  +  [-(jy  +  U>xU}z]Hr  +  [-(^y  +  f^lVr'^Zg  (A.17) 

The  net  torque  generated  is  given  by: 

TL^fl  +  f^  (A.18) 


In  the  derivation  of  the  torque  generated,  the  mass  centers  were  allowed  to  be  different 
distances  from  the  origin.  The  device  was  not  symmetric  about  the  ajy,  j/y,  or,  Zg  axes.  With 
modern  fabrication  facilities  and  the  ability  to  laser  tune  devices,  the  mass  centers  may  be 
assumed  to  be  located  an  equal  distance  from  the  origin,  resulting  in  symmetry  about  the 
sensitive  axis.  In  addition,  the  mass  of  the  two  members  will  be  the  same.  Under  these 
assumptions,  the  following  simplifications  may  be  made: 


II 

II 

(A.19) 

II 

II 

(A.20) 

bi  —  br  —  b 

(A.21) 
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rriL  =  nir  =  m 


(A.22) 


The  net  torque  ff„t  is: 

Ttot  =  2m[[(d;j;  +  UyU}z)b‘^  +  +  ^z)Hb+  {Uyio^  —  u:^)H‘^]Xg 

+[(r-  +  9,sin{ut))[-[2Ljyil  +  [uy  +  Lj^u}^]Q.]sin{u}t)  -  2Q.UyUcos{ut) 

+  niyjr]  +  [u^ujy  -uj:,]Hh]yg 
+[(r  +  nsin(a;t))[-[2cj^l7  +  {u>xt^y  +  u>^)^i\sin{ut)  —  2^u)  ^u}cos(ojt)] 

-[LJ^LOy  -  Lj^]b‘^]zg]  (A.23) 

The  above  case  allows  for  a  displacement  of  the  mass  center  of  the  system  in  the  yg 
direction.  This  equation  may  be  further  simplified  by  assuming  a  uniform  mass  distribution 
in  the  vibrating  members.  The  result  of  this  assumption  is  a  zero  displacement  of  the 
position  vectors  in  the  yg  direction  and  symmetry  about  a  line  joining  the  two  mass  centers. 
To  account  for  this  in  the  equations,  let  6  =  0.  The  output  torque  then  becomes: 

ff„t  =  2m[[{uyU;^-u^)H^]xg 
+[(r  +  ^sin{u!t))[—[2UyCl  +  [Cjy  +  u};iU>x\^]sin{uJt) 

-2UuyUcos{ujt)  -  +  Uy]r]yg 

+[(r  +  Qsin{u>t))[—[2u>^il  +  +  d;^)fi]sm(a;t) 

— 2flu;^u;cos(a;t)  —  [u^Uy  +  0Jz]'f‘]zg]  (A.24) 


A-6 


This  is  the  final  form  for  the  output  torque  that  will  need  to  be  implemented  for  arbitrary 


rotations.  In  order  to  verify  the  single  axis  rotation  equation  presented  in  Chapter  3,  let 
the  input  rotation  be  only  about  the  Zg  axis,  and  let  there  be  no  change  in  the  maximum 
amplitude  of  the  sustained  oscillation,  so; 


=  ca,,  =  =  0 


(A.25) 


the  output  torque  becomes: 


ffo-t  =  2m[(r  +  Q.sin{u)t))[—u);i[r  +  Qsin{ojt))  —  2^lu^ojcos{ujt)]]zg  (A. 26) 


This  is  identical  to  the  result  from  Chapter  3  for  an  input  rotation  rate  of  Ug. 
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Appendix  B.  Planar  Gyroscope  Output  Derivations 


The  full  derivation  of  the  planar  gyroscope  torque  equations  is  presented  for  an 
arbitrary  input  rotation.  The  position  vectors  to  the  mass  centers  are  shown  in  Figure 
B.l.  The  four  position  vectors  Ri,R2,R3.,  and  R4  are  defined  by  Equations  (B.1-B.4). 


Figure  B.l  Block  diagram  for  generic  torque  derivations  planar  gyroscope. 


Ri  = 


D 

~R 


0 


(B.l) 
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(B.2) 


(B.3) 


(B.4) 


The  sustained  oscillation  is  defined  by  Equation  (B.5): 


Ay  —  Q,sin{u>t)yg 


(B.5) 


Let  the  input  rotation  be  defined  by  Equation  (B.6): 
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This  results  in  a  net  rotation  rate  with  respect  to  inertial  space  given  in  Equation  (B.7): 


ujy  +  uQ.cos{u}t) 


(B.7) 


The  inertial  accelerations  of  each  of  the  four  position  vectors  are  then  calculated  and  shown 
in  Equations  (B.8-B.11): 


Ai  =  [RlOz  +  {u>y  +  ujQcos{ojt))(— RiJx  —  Dojy  —  DuiQcos^ojt))  +  Du'^]xg 
+[-D6j^  -  RujI  —  DuxOJy  —  Du>j;uQ,cos{u;t)  —  Ruj^]yg 
+[—Rojx  —  Dojy  -  D6jQ,cos{ut)  +  Du'^Q.sin(ut)  —  —  Roj^Uy 

—Roj2U>Q,cos(ujt)]zg  (B.8) 

A2  =  [—Ru^z  +  (wy  +  u>Qcos{u)t)){RoJz;  —  Duy  —  DuQcos^ojt))  +  D(jj^]xg 
+[— +  R(x)l  —  Du^Uy  —  Duz:U^cos{ujt)  +  Rw1]yg 
+[i2(hj;  —  Duy  —  DuQ,cos{u;t)  +  DuP^sin{ut)  —  +  RuJzOJy 

-\-RLJzUj^cos(ujt)]zg  (B.9) 

A3  =  [—RUz  +  i'^y  +  U!^COs{(jtjt)){RUx  +  DUy  +  DujUcOs{ujt))  —  Du)^]Xg 
+  [Du;z  +  Rl^l  +  DUz^iVy  +  DUz:ljJ^COs(jjjt)  +  RLjl]yg 
+[RLbz:  +  Duy  +  DojQcos{ivt)  —  Duz^Q,sin(tjjt)  +  Doj^^Jz  +  RuJz^y 

■]-RijJzijjQ,cos{(jji)\zg  (B.IO) 

A4  =  [RuJz  +  (wy  +  i^^cos{LJt)){—Rux  +  Dojy  +  Du;^cos{u)t))  —  Du^]xg 
+  [Du}z  —  Ri^l  +  Du)z;iZ)y  +  DUz:UQ,COs(ut)  —  Rijl]yg 
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+[— +  Dujy  +  DCj£Lcos{(jjt)  —  Du'^Q,sin{LJt)  +  Doj^ui^  —  Ruj^Wy 

—Rbj^ujQ,cos{ujt)\zg  (B.ll) 

Using  the  derived  accelerations  and  the  definition  of  torque,  the  torques  generated  by  each 
of  the  mass  centers  are  calculated  and  given  in  Equations  (B.12-B.15): 

Ti  —  +  DRjhy  +  DRu)Q,cos(ojt)  -  DRuj'^9,sin(Ljt) 

+DRll!xU>;,  +  R'^ojyio^  +  R?(jj!,u>^cos(ut)]x  g 

+[—RDu}x  —  —  D^(hQ.cos{cot)  +  D'^u)^Qsin{ut) 

—  DRu>yU2  —  DRu)^uQ,cos{ujt)]yg 
+[-D^w^  -  DRujI  -  D^UxUy  -  D^u)xOjQcos{ut)  +  R^oj^ 

+(u)y  +  u^cos{ut){- R?ujx  —  DRwy  -  DRu>Q,cos{u}t))]zg]  (B.12) 

T2  =  m[[i2^Wj.  —  DRCjy  —  DRu>^cos{ojt)  +  DRu>^^lsin{ut) 

—DRio^uj^  +  R^UyUJz  +  i2^a;^a;Ocos(wt)]£j 
+[RDCjg;  —  D^ojy  —  ujQ,cos{(jjt)  +  D'^uj^ilsin{u)t) 

—D^u>x(jJz  +  DRWyOJ^  +  DRuj2W^cos{u!ty\yg 
+  DRojI  —  D^ojj;Uy  —  D'^LJx^Qcos{u)t)  +  R^Ulz 
+(u;j,  +  ojncos{ojt){—R‘^Ux  +  DRuy  +  DRu^cos{u>t))]zg]  (B.13) 

T3  =  'm[[R'^Ux  +  DRuy  +  D  RxvQ,cos(u>t)  —  DRu'^9,sin{ujt) 

+  DRUxOJ;z  +  R?bOyOJ2  +  R?U}zbjO.COs{ujt)\Xg 
-{^[—RDujx  —  —  D^u^cos{ut)  +  D^u^ilsin{ujt) 
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—D^UxOJz  —  DRujyU^z  —  DRuj^,uQ,cos[ujty\yg 
-\-[—D'^u^  —  DRuil  —  D^ujxOJy  —  D'^Ux^^c.os{u)t)  +  B?Cj^ 

+(a;j,  +  ujQ.cos{ut){- R^lJx  —  DRuy  —  DRu)^cos(LL>t))]zg]  (B.14) 

T4  =  m[[R^ui^  —  DRuy  —  DRujilcos{ut)  +  DRjjj‘^Q,sin{ujt) 

—DRu^u^  +  R^ojyU)^  +  R'^(jj^LjQ,cos{u:t)]xg 

+[RDd!x  —  —  D^u}Q,cos(u>t)  +  D'^Lj'^£lsin{u>t) 

—D'^UxOJz  +  DRjbjyU^  +  DRu;:Uj^cos{u}t)\yg 

+[— +  DRul  —  D'^ujxOJy  —  D‘^UxU}Q,cos{ut)  +  R?u>^ 

+{ojy  +  (jj9.cos{u!t){—R'^oj^  +  DRojy  +  DRuQ,cos{u>t))]zg]  (B.15) 

The  total  torque  is  the  sum  of  the  torques  generated  from  each  of  the  mass  centers.  The 
resulting  net  torque  is  given  by  Equation(B.16): 

T  =  47n[[R^{u)x  +  i^yijJz  +  ojz<^^cos{ijt))]x  g 
+[D'^{—Uy  —  —  djQcos{ut)  +  uj'^ilsin{uzt))]yg 

+  [cj^(i?^  —  D^)  —  D'^ljJj:(Vy  —  D^iX)z;UQ,COs{Ljt) 

—  R?bJz:U)y  —  R?U>3;UzUcOs{ujt)]Zg  (B.16) 

The  net  torque  simplifies  to  Equation  (B.17);  this  is  identical  to  the  result  derived  in 
Chapter  3. 


T  =  Am[[R^CjJzU}ilcos{u;t)]x  g 
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-\-[—D'^u:D,cos{u>t)  +  D'^uy‘'Q,sin{u>t))\yg 


(B.17) 
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Appendix  C.  Ksys  Matlab  Script  File 


E=160e6; 

rho=2330; 

num_brace=4 ; 

wf inger=2e-6; 

scale=.3714; 

lplate=50e-6; 

wplate=100e-6; 

sigr=input( ’residual  stress  of  layer  :  ’); 
sigr=sigr*le6 ; 

thick=input( 'Layer  Thickness  in  microns 
thick=thick*le-6 ; 

wsup=input(’ Support  Width  in  microns  , 
wsup=wsup*le-6; 

lsup=input( ’Support  Length  in  microns  :’); 
lsup=lsup*le-6 ; 

If ing=input( ’Finger  Length  in  microns  :’); 

If ing=lf ing*le-6 ; 

numfingers=input(’ Number  of  Fingers  :’); 
kl=(wsup*thick) *sigr/lsup ; 
k2=(E*thick*(wsup)''3)/(lsup“3) ; 
ksys=num_brace*(kl+k2) ; 

Mp=wplate*lplate*thick*rho ; 

Mf =numf ingers* (If ing*thick*wf inger) ; 

M=Hp+Mf ; 

Msup=wsup*thick*lsup*rho ; 

Fr= (1/ (2*pi) ) * (ksys/ (M+scale* (Msup) ) ) “ . 5 ; 
f ile=fopen( ’freq’ , ’a’) ; 

fprintf (file, ’Output  for  Tunning  Fork  Gyroscope  \n’); 
fprintf  (file, ’Youngs  modulus:  '/,10.3e  \n’ ,  E)  ; 
fprintf  (file, ’Rho  used:  '/,10.3e  \n’ ,  rho)  ; 
fprintf  (file, ’number  of  braces:  '/,10.3e  \n’ ,  num_brace)  ; 
fprintf  (file, ’finger  length:  '/,10.3e  \n’ ,  If  ing)  ; 
fprintf  (file, ’Support  length:  '/.lO.Se  \n’ ,  Isup)  ; 
fprintf  (file, ’Support  width:  y,10.3e  \n’ ,  wsup) ; 
fprintf  (file, ’Layer  Thickness:  */,10.3e  \n’,  thick); 

fprintf  (file, ’scale  factor  for  spring  constants:  */,10.4e  \n’ ,  scale); 
fprintf (file, ’kl  k2  ksys  Mp  Mf 

M  Msup  Fr  \n’); 

fprintf  (file,’ y.S.Se  y,5.3e  y.5.3e  y.5.3e  y.5.3e  y.5.3e  y.8.3e  y,5.3e\n\n\n’ 
,  kl, k2, ksys, Mp,Mf,M, Msup, Fr) ; 
fclose(’all’) ; 
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